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ABSTRACT 

We present the results from an X-ray and optical study of a new sample of eight ex- 
treme luminosity ultraluminous X-ray source (ULX) candidates, which were selected 
as the brightest ULXs (with Lx > 5 x 10 40 erg s^ 1 ) located within 100 Mpc identified 
in a cross correlation of the 2XMM-DR1 and RC3 catalogues. These objects are so 
luminous that they are difficult to describe with current models of super-Eddington 
accretion onto all but the most massive stellar remnants; hence they are amongst the 
most plausible candidates to host larger, intermediate- mass black holes (IMBHs). Two 
objects are luminous enough in at least one observation to be classed as hyperluminous 
X-ray source (HLX) candidates, including one persistent HLX in an SO galaxy that 
(at 3 x 10 41 erg s _1 ) is the second most luminous HLX yet detected. The remaining 
seven sources are located in spiral galaxies, and several appear to be closely associated 
with regions of star formation as is common for many less luminous ULXs. However, 
the X-ray characteristics of these extreme ULXs appear to diverge from the less lu- 
minous objects. They are typically harder, possessing absorbed power-law continuum 
spectra with L ~ 1.7, and are potentially more variable on short timescales, with data 
consistent with ~ 10-20 per cent rms variability on timescales of 0.2-2 ks (albeit at 
low to moderate significance in many datasets). These properties appear consistent 
with the sub-Eddington hard state, which given the observed luminosities of these 
objects suggests the presence of IMBHs with masses in the range 10 3 — 10 4 M Q . As 
such, this strengthens the case for these brightest ULXs as good candidates for the 
eventual conclusive detection of the highly elusive IMBHs in the present-day Universe. 
However, we caution that a combination of the highest plausible super-Eddington ac- 
cretion rates and the largest permitted stellar black hole remnants cannot be ruled 
out without future, improved observations. 
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1 INTRODUCTION 

Although the existence of a population of black holes 
with masses between those in Galactic black hole X- 
ray binaries (BHXRBs; ~ 10 Mq) and in active galac- 
tic nuclei (AGNs; > 10 5 Mq) has long been postu- 
lated, these intermediate-mass black holes (IMBHs) have 
remained observationally elusive, with dynamical evidence 
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for such objects in large glo bular clusters still th e sub- 
ject of some dispute (e. g. Gebhardt. Rich fc Hoi 120051 ; 



l Andersonfcvander M^ T A Yet such objec t"- 
pear necessary to explain the initial seeding of the su- 
permassive bla ck holes (SMB Hs) now present in most 
galaxy nuclei (|Volonteril I2010T ). whose mass was subse- 
quen tly built up either through the accretion of gas 
(e.g. lYu fc Tremaine 20021). or through hierarchical black 
hole mergers (e.g. ISchneider et al.ll2002l ). It is likely that 
the processes forming these seeds - for example, the col- 
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lapse of Population III stars, or dynamical interactions in 
dense stellar environments - would have l eft relic IMBHs 
in the halos of present-day galaxies (e.g. | Madau fc Reesl 
l200ll ; lEbisuzaki et alj l200ll : llslam et all 120041 ). Indeed, the 
latter process may still be ongoing in som e young, dense 
stellar clusters (jPortegies Zwart et al.l 12004 ). adding to the 
primordial IMBH population. Furthermore, the tidal strip- 
ping of merging satellite galaxies containing nuclear black 
hole s may also add to this population in large galaxy haloes 
fe.g. iBellovarv et alj|20ich . It is therefore puzzling that the 
existence of such IMBHs remains to be confirmed. 

The detection of a population of unusually luminous 
X-ray sources - in excess of 10 39 erg s _1 - located out- 
side the nuclei of nearby galaxies initially led to sugges- 
tions that these objects could harbour the long sought- 
after IMBHs (e.g. lColbert fc Mushotzkvlll999l ). Early Chan- 
dra & XMM-Newton spectral re sults for these ultr alumi- 
nous X-ray sources (ULXs; see iFeng fc Sorial 120111 for a 
recent review) appeared to support this diagnosis with 
the detection of a pparent cool disc signatures, indica- 
tive of IMBHs fe.g. iKaaret et al.l 120031 ; IMiller et all 120031 ; 
iMiller. Fabian fc Miller] 12004). However, subsequent analy- 
ses of high quality XMM-Newton data have instead revealed 
spectral and temporal characteristics in many ULXs that ap- 
pear inconsistent with the sub-Eddington accretion regime 
expected for IMBHs at typical ULX luminosities. Instead, 
these analyses support models where ULXs are powered by 
super-Eddingto n accretion on t o stellar remnant black holes 
(< 100 Ma , iRobertsI 120071 : Idadstone. Roberts fc Done 



20111 ). 



tsl 120071: 

Si 



2009; iMiddleton et a.1.1 I2OIII ; iMiddleton. Sutton fc Roberfc 



Further observational support for ULXs as stellar-mass 
objects comes from the close spatial and temporal as- 
sociation of ULXs with ongoing star for mation, as most 
notably seen in the Cartwhe el galaxy ([Gap et al.l 12003 : 
lKindl2004 iMapelli et"ai"1l200Sl ). suggesting they are a type 
of high-mass X-ray binary (HMXB) which naturally pro- 
vides them with the la rge mass reservoir required to sus- 
tain ULX luminosities (jRappaport. Podsiadlowski fc Pfahll 
120051 ). It also comes from the power- law form of the 
X-ray luminosity function (XLF) for HMXBs in nearby 
star forming galaxies, that is smooth across the transi- 
tion from standard HMXBs to ULXs, but then breaks 
at ~ 2 x 10*° erg s" 1 (iGrimm Gilfanov fc Sunvaevl 12003" ; 
ISwartz et alll2004l ; ISwartz et al.ll201ll ). This break is impor- 
tant, as its luminosity is only ~ 10 per cent of the Edding- 
ton luminosity for a ~ 10 3 Mq IMBH; a similar turn-off at 
such a low Eddington fraction is not seen in other accretion- 
powered populations, making it difficult to reconcile with a 
population of IMBHs, but is again supportive of most ULXs 
as a population of super-Eddington BHXRBs. 

However, there remain some plausible IMBH candi- 
dates in the ULX population. Recent work bv lSwartz et al.1 
l|201ll ) has demonstrated that an extrapolation of the local 
XLF out to 100 Mpc underpredicts the observed numbers 
of high luminosity ULXs, particularly those that possess lu- 
minosities greater than 10 41 erg s _1 , suggesting that such 
objects originate in a separate physical population to the 
HMXBs that likely underlie most ULXs. Indeed, while a 
combination of massive stellar black h oles, possibly formed 
in low metallicity regions (up to 80 Mm; IZampieri fc Roberts! 
2009. Belczvnski et alJl201CT ) and close to maximal radiation 



from super-Eddington accretion (~ 20 Z/Edd from a face-on 
disc; lOhsuga fc Mineshigell201ll ) could explain ULXs up to 
~ 10 41 erg s _1 , but objects exceeding this luminosity appear 
to require a different physical scenario. It is these hyperlu- 
minous X-ray sources (HLXs; Lx > 10 41 erg s _1 ) that may 
present us with the best evidence for IMBHs to date. 

Very few bona fide HLXs are known: the b est candidates 
Matsumoto et al.ll200ll). ESQ 243-49 



include M82 X-l (e.j_, .. ... _. 

HLX -1 (|Farrell et al.ll2009l ), Cartwheel N10 (iPizzolato et al.1 



120101 ) and CXO J122518. 6+ 144545 (|jonker et al.l l201Ch . 
Most observations have focussed on the nearest HLX, M82 
X-l at d ~ 4 Mpc, and the most luminous, ESO 243-49 
HLX-1 at Lx.poak ~ 10 42 erg s" 1 . M82 X-l is notable for 
several reasons, inclu ding its co-location with a massive, 
young stellar cluster (|Portegies Zwart et al.ll2004l ). the de- 
tection of a 50 - 100 mHz QP O in its power spectrum 
l|Strohmaver fc Mushotzkvl 120031), a reported 62-day peri- 
odicity in its X-ray light curve (jKaaret fc Fend 120071 ) and 
plausible transiti ons between a hard state and a thermal 
dominant state (|Feng fc Kaaretl |2010j), all of which con- 
tribute to the picture of this object as an excellent IMBH 
candidate. ESO 243-49 HLX-1 is most remarkable for its 
high peak luminosity and regular outbursts, which appear 
to possess a fast rise/expon ential decay profile and to repeat 
on timescales of ~ 1 year (JLasota et al.ll201ll ). This object 
also appears to mimic BHXRB behaviour in its hardness- 
intensity diagram as it transits throu gh its outbursts, indic a- 
tive of sub-Eddington state changes (|Servillat et al.ll201lT ). 

It is evident that the most luminous ULXs are rare, 
and little is known about them as a class; yet we have sum- 
marised good arguments as to why these are potentially im- 
portant targets in the ongoing search for IMBHs. One way 
to learn more of this population is to extend the number 
of known objects using archival searches. Here, we report 
a study of te n candidate extreme ULXs selected from the 
IWalton et al.l (2011b) catalogue of ULXs detected in the 
2XMM-DR1 survey, that display observed luminosities well 
in excess of the break in the XLF, including two new HLX^J. 

The paper is laid out as follows. In section 2 we de- 
scribe the selection of the new sample, and the reduction 
of the recent optical and X-ray data associated with them. 
Section 3 details the imaging, spectral and timing results 
from all available XMM-Newton and Chandra datasets for 
the sources, and describes an optical follow up study. The 
characteristics of these objects are then compared to those 
of the less luminous bulk of the ULX population in section 
4, before the implications of the findings are discussed with 
respect to the nature of the objects in section 5, and the 
paper is concluded in section 6. 



2 SAMPLE SELECTION & DATA REDUCTION 

2.1 Sample selection 

IWalton et al.l (2011b) present a catalogue of 470 candi- 
date ULXs in 238 galaxies, detected in archival XMM- 
Newton observations, produced from a cross-corr elation of 
the RC3 catalogue of galaxies (jde Vaucouleurs et al.lll991i ) 



1 This work builds and improves on th e early results from this 
sample, published in Su tton et al.1 I I2011I) . 
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and the 2XMM-DR1 catalogue <| Watson et all 120091 ). This 
was reduced to a small sample of the most luminous, rel- 
atively nearby ULXs by filtering on luminosity (Lx > 
5 x f0 40 erg s _1 ) and distance (d < fOO Mpc). A sample 
of 12 sources were found to match these criteria. Additional 
quality checks of the XMM-Newton data for the selected 
sources were then performed, and 2 of the 12 sources were 
excluded from the final sample. We excluded a ULX candi- 
date detected near to the cD galaxy NGC 4889, in the galaxy 
cluster Abell 1656, where we could not separate sufficient 
source counts from the strong diffuse emission in which the 
candidate ULX was embedded to permit detailed study. The 
other excluded source was M82 X- l, as it has been exten- 
sively studied by other authors (e.g . Matsu moto et al.ll200ll ; 
iKaaret. Simet fc Lang|l2006l : iFeng fc Kaaredl2010h . and the 
XMM-Newton data is somewhat complicated by the pres- 
ence of a strong ISM component around the position of the 
ULX, and bright and variable nearby object^ 2 ]. 

The 10 remaining objects are listed in Table (TJ along 
with some characteristics of their presumed host galaxies. 
The X-ray sources are assumed to be at the distance of 
the host galaxy, and for the majority of sources cosmology- 
corrected distances are used (Hp = 71 km s _1 Mpc -1 , 
Q M = 0.27, Q A = 0.73. iKomatsu et aT|l2009h , calculated 
using redshifts from lde Vaucouleurs et al.l 1 1991T ). An excep- 
tion to this is Src. 8, whose host galaxy is sufficiently nearby 
that its peculiar motion could dominate its measured reces- 
sion velocity , so we instead use the distance of 14.9 Mpc from 
iTullvl (|l988T ). Table Q] also provides a luminosity for each ob- 
ject, calculated using this distance and the 2XMM 0.2-12 
keV flux. It is immediately evident that the extreme lumi- 
nosity sample presented here includes 3 candidate HLXs - 
Src. 1, Src. 4 and Src. 7 (add itionally, Src. 3 had previo usly 
been identified as a HLX bv lDavis fc Mushotzky||2004l . see 
section I3.1.1J I. 

In addition to the data sets identified in the 2XMM 
catalogue, the XMM-Newton science archivqj and Chan- 
dra data archivqjwere searched for further observations cov- 
ering the objects. All the archived XMM-Newton EPIC and 
Chandra ACIS observations potentially covering the posi- 
tions of the 10 candidate ULXs were extracted from the 
archive, and examined for further detections of the 10 ob- 
jects. Details of all relevant observations are shown in Ta- 
ble [2] with the exceptions of the XMM-Newton observations 
0300530401 and 0301651201 (plus several additional short 
observations of Src. 6, see caption for details). These were 
pointed in the proximity of Src. 6 and Src. 10 respectively, 
but were not used as they contained no useful data due to 
the presence of heavy background flaring throughout the 
duration of the EPIC exposures. In all cases, the ULX in 
question was detected. All but one of the sources (Src. 7) are 
present in multiple archived observations, including six with 
archival Chandra data. We were awarded additional Chan- 
dra time during cycle 12 to observe the other four sources 
(proposal ID 12620389, PI Sutton; results for Src. 2, Src. 



4 and Src. 7 have been obtained and are included here, but 
an observation of Src. 8 awaits scheduling at the time of 
writing). 



2.2 XMM-Newton data reduction 

The XMM-Newton data were reduced and data products 
were extracted using the standard tools available in the 
XMM-Newton SAtQ (version 9.0.0). Severe background flar- 
ing resulted in multiple exposures being obtained in a num- 
ber of the observation^. However, in all cases only one of 
the exposures contained a significant length of good time. 
In each case we identified the sole exposure containing use- 
ful science data per observation and focused our analysis 
upon it. In two observations of Src. 6 the solar proton flar- 
ing was severe enough to render the data from whole instru- 
ments unusable, in particular the PN detector in observation 
0300530301 and the MOS1 detector in 0300530601. How- 
ever, background flaring was present to some significant level 
in almost all the XMM-Newton observations, with the only 
exceptions being observations 0300530301 and 0112980101 
of the ULX candidates Src. 6 and Src. 9 respectively. Good 
Time Interval (GTI) files were therefore created to allow 
high background epochs to be filtered out of the data, based 
on high energy (10-15 keV) full field light curves. The limit- 
ing count rate used to define the GTIs varied between obser- 
vations, dependent on the best balance between minimising 
the source data loss whilst maximising the excision of back- 
ground flares, with typical thresholds of 1.0-1.5 cts s _1 for 
PN fields and ~ 0.6 cts s" 1 for MOS fields. 

Source spectra and light curves were then extracted 
from circular apertures centred on the ULX candidate. 
Source extraction region radii were typically between ~ 15- 
40 arcseconds (which corresponds to an on axis fractional 
encircled energy of ~ 70-90 per cent at 1.5 keV), in order to 
maximise source counts whilst avoiding contamination from 
nearby sources. In three cases the source was positioned close 
to a chip gap, which would have severely restricted the ra- 
dius of a circular aperture. In each of these cases elliptical 
source regions (with semi-major axes equivalent to the cir- 
cular aperture radii) were used to extract the source data. In 
both the MOS1 data from Src. 3 in the 0022340201 dataset, 
and the PN data for Src. 5 in observation 0147610101, the 
ellipse was aligned with the nearby chip gap; however in the 
PN data for Src. 2 in 0093630101, the source was sufficiently 
far off-axis that the ellipse was aligned with the direction of 
maximum point spread function (PSF) extent. Larger circu- 
lar apertures were used to provide source background data 
in all cases. These were positioned on the same chip, ap- 
proximately the same distance from the readout nodes if 
possible, and were sized to include a minimum of ~ 50 per 
cent more area on the detector than source apertures. In 
cases where nearby sources placed limitations on the size of 
a background region, multiple circular apertures were used. 



2 Of the other notable HLXs mentioned in the introduction, 
ESO 243-49 HLX-1 was not included in the sample as its host 
is not in the RC3 catalogue, and both Cartwheel N10 and CXO 
J122518.6+144545 lie beyond the 100 Mpc limit of this sample. 

3 http://xmm.esac.esa.int/xsa/ 

4 http://cxc.harvard.edu/cda/ 



5 http://xmm.esac.esa.int/sas/ 

6 This affected both MOS and PN detectors during observation 
0145190201 of Src. 8, and MOS detectors alone were affected dur- 
ing observations 0200780101, 0147610101 and 0145190101 of Src. 
1, Src. 5 and Src. 8 respectively. 
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Table 1. The extremely luminous candidate ULX sample. 



Source ID 



2XMM source 



Host galaxy Galaxy type a Distance 6 Separation 



N H d 



ix,u 



Src. 
Src. 
Src. 
Src. 
Src. 
Src. 6 
Src. 7 
Src. 8 
Src. 9 
Src. 10 



2XMM 
2XMM 
2XMM 
2XMM 
2XMM 
2XMM 
2XMM 
2XMM 
2XMM 
2XMM 



J011942. 
J024025. 
J072647. 
J120405. 
J121856. 
J125939. 
J134404. 
J151558. 
J163614. 
J230457. 



7+032421 
6-082428 
9+854550 
8+201345 
1+142419 
8+275718 
1-271410 
6+561810 
0+661410 
6+122028 



NGC 470 
NGC 1042 
NGC 2276 
NGC 4065 
NGC 4254 
NGC 4874 

IC 4320 

NGC 5907^ 

MCG 11-20-19 

NGC 7479 



SA(rs)b 

SAB(rs)cd 

SAB(rs)c 

E 

SA(s)c 

cDO 

SO? 

SA(s)c: edge-on 

Sa 

SB(s)c 



32.7 
18.9 
33.3 
87.9 
33.2 
99.8 
95.1 
14.9 
96.2 
32.8 



33 


3.09 


15.3 + 0.8 


96 


2.61 


5.1 + 0.2 


15 


5.52 


6.4 + 0.3 


18 


2.40 


13.1 + 1.7 


103 


2.70 


6.0 + 0.3 


57 


0.89 


5.8 + 1.3 


18 


5.04 


28.2 + 2.4 


102 


1.38 


5.6 + 0.1 


16 


3.14 


6.5 + 1.2 


68 


5.08 


7.1 + 0.3 



Notes: a Galaxy mor phology from |de Vaucpuleu rs et al.l (1199 lh . b Host galaxy distances in Mpc. These are cosmo logy-c o rrecte d distances 
using redshifts from Ide Vaucouleurs et al.l (J1991J), with the exception of NGC 5907 where the distance from iTullvl 1(1988') was used. 
c Separation of t he candidate ULX positio n from the nucleus of its host galaxy, in units of arcseconds. d Galactic column densities 
interpolated from lDickev fc Loc kman (1990J), in units of 10 20 cm -2 . e Maximum detected 0.2-12 keV luminosity of the ULX candidate, 
based on the 2XMM broad band fluxes and the quoted distance, in units of 10 40 erg s~ . * The part of NGC 5907 west of its dust lane, 
that hosts the ULX candidate, is also known as NGC 5906. 



The final data products used in the analysis (i.e., spec- 
tra and light curves) were then extracted from the source 
apertures, using the GTI file and standard event patterns 
for the detectors (PATTERN < 4 for the PN, PATTERN < 12 
for the MOS). Spectral extractions also used the conserva- 
tive FLAG = filter, although we alleviated this somewhat 
for the light curve extractions by using the less restrictive 
#xmmea_ep or #xmmea_em filters for the PN or MOS detectors 
respectively instead. The appropriate response matrices for 
the spectral analysis were automatically generated during 
the spectral extraction in the SAS task xmmselect. 



2.3 Chandra data reduction 

Chandra data were reduced using ClAqj version 4.1.1, and 
calibration files from the Chandra CALDB version 4.1.0. Cir- 
cular apertures with radii of 4 arcseconds (corresponding to 
an on axis fractional encircled energy of > 99 per cent at 
1 keV) were used to extract the source data, with the only 
exception being in the case of Src. 3, as this ULX candi- 
date was resolved by Chandra into multiple objects sepa- 
rated by less than the typical extraction region radius (see 
section [37T7TJ . Instead, apertures of ~ 1 arcsecond radius 
(~ 90 per cent on axis fractional encircled energy at 1 keV) 
were used to avoid overlapping the source regions. Annu- 
lar apertures, centred on the source, but with inner radii 
outside of the source region were used to characterise the 
background for each source. In the case of Src. 3, a sin- 
gle background annulus surrounding the resolved sources 
was used to avoid cross-contamination of the background 
regions by the neighbouring point sources. The data prod- 
ucts were then extracted from the archived level 2 event files, 
to include all standard good event detections with energies 
between 0.3-10 keV, within the source data extraction aper- 
tures. The Chandra spectra were extracted using the CIAO 
script SPECEXTRACT, which also produces the appropriate 
response matrices for spectral analysis, and the light curves 
were extracted using dmextract. 



2.4 Optical data 

This work also utilises optical data to further probe these 
extreme objects. As an initial step we retrieved Digitised 
Sky Survey version 2 (hereafter DSS2) data from the ESO 
archivqfl, within a 5 arcminute square centred on the host 
galaxy of each ULX candidate. DSS2-blue images were used 
preferentially when available, but due to its limited coverage 
(only 45 per cent of the sky), the use of the alternative all- 
sky DSS2-red images was necessary for Src. 2 and Src. 7. 

An examination of the Hubble Legacy Archive (HLA jj 
revealed that high resolution space-based optical imaging 
from HST was available for four of our objects. Enhanced 
HLA data products covering the source region were available 
for Src. 3, Src. 6 and Src. 8, and two data sets processed us- 
ing the standard HST single image calibration pipeline were 
available for Src. 10. All this data was downloaded from 
the archive for further study; we tabulate the basic obser- 
vation details in Table [3] As an initial step the astrometry 
of the images was improved using the iraf tools CCFIND, 
CCMAP and CCSETWCS and objects common to the HST im- 
ages and either the 2MASS (Src. 8) or USNO (Src. 3, Src. 
6 and Src. 10) catalogues. Where individual counterparts to 
the X-ray sources were detected, aperture photometry was 
performed using GAIA and the appropriate Vega magnitude 
zero points from either ISirianni et alj (|2005l ) for ACS/WFC 
observations, or the WFPC2 Data Handbook 10 !. Corrections 
to a standard 0.5 arcsecond aperture were calculated using a 
mean magnitude correction from ten bright point sources in 
each field, and these were extrapo lated to an infinite r adius 
aperture, again using values from ISirianni et al.l (|2005l ) and 
the WFPC2 Data Handbook. 

In addition to the study of archival imaging data for 
the sources, a low dispersion Gemini CMOS spectrum was 
obtained to further constrain the nature of an optical coun- 
terpart to Src. 6. The data were reduced using standard 
tools in the gmos sub-package from the Gemini package in 



7 http://cxc.harvard.edu/ciao/ 



8 http://archive.eso.org/dss/dss 

9 http://hla.stsci.edu/ 

10 http : //www . stsci . edu/hst/wf pc2/documents/handbook/WFPC2_DHB . html 
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Table 2. Observation details 



Source ID 


Observatory 


ObsID Q 


Date 


t b 


Total counts c 


Off-axis angle d 










(yyyy/mm/dd) 


(ks) 


(0.3-10.0 keV) 


(arcmin) 


Src 


. 1 


XMM-Newton 


0200780101 


2004/01/24 


4.0/3.8/3.6 


850 ± 30 


5.25 






Chandra 


10562 


2008/11/12 


28.1 


610 ±20 


0.29 






XMM-Newton 


0601670101 


2009/06/27 


40.1/53.9/47.1 


2550 ± 50 


1.10 


Src 


. 2 


XMM-Newton 


0093630101 


2001/08/15 


11.6/15.0/ - 


920 ± 30 


14.41 






XMM-Newton 


0306230101 


2006/01/12 


- /53. 1/52.1 


2860 ± 60 


12.18 






XMM-Newton 


0553300401 


2009/08/12 


42.2/48.4/ - 


4310 ± 70 


14.38 






Chandra 


12988 


2010/10/09 


29.0 


2770 ± 50 


0.32 


Src 


. 3 


XMM-Newton 


0022340201 


2001/03/16 


- /52.3/52.3 


2070 ± 50 


6.08 


Src. 


3a e 


Chandra 


4968 


2004/06/23 


45.6 


320 ± 20 


1.13 


Src. 


3b e 


Chandra 


4968 


2004/06/23 


45.6 


190 ± 10 


1.10 


Src. 


3c e 


Chandra 


4968 


2004/06/23 


45.6 


120 ± 10 


1.06 


Src 


. 4 


XMM-Newton 


0112270601 


2003/01/02 


- /6. 5/6.5 


75 ±9 


2.10 






XMM-Newton 


0112271001 


2003/06/09 


- /6.9/6.9 


140 ± 10 


0.48 






XMM-Newton 


0112271101 


2003/06/30 


3.1/3.6/3.6 


210 ±10 


0.60 






Chandra 


12990 


2011/07/13 


5.0 


57 ±8 


0.30 


Src 


. 5 


XMM-Newton 


0147610101 


2003/06/29 


13.1/17.3/19.4 


1770 ± 40 


2.75 






Chandra 


7863 


2007/11/21 


5.1 


270 ± 20 


1.84 


Src. 


6/ 


Chandra 


556 


1999/11/04 


9.6 


60 ±10 


2.28 






XMM-Newton 


0124711401 


2000/05/29 


11.0/16.8/16.8 


350 ± 80 


0.42 






XMM-Newton 


0153750101 


2001/12/04 


17.0/21.0/21.0 


580 ± 90 


2.69 






XMM-Newton 


0300530701 


2005/06/06 


16.0/23.6/24.3 


430 ± 90 


2.14 






XMM-Newton 


0300530601 


2005/06/07 


14.3/ - /22.7 


410 ± 80 


2.21 






XMM-Newton 


0300530501 


2005/06/08 


15.5/24.8/24.9 


480 ± 90 


2.31 






XMM-Newton 


0300530301 


2005/06/11 


- /30.4/30.4 


230 ± 70 


1.60 






XMM-Newton 


0300530101 


2005/06/18 


15.7/18.0/18.0 


480 ± 90 


1.42 






Chandra 


9714 


2008/03/20 


29.7 


40 ±10 


1.86 






Chandra 


10672 


2009/03/15 


28.5 


110 ±10 


1.93 


Src 


. 7 


XMM-Newton 


0092970801 


2003/01/04 


3.5/5.3/5.2 


390 ± 20 


1.10 






Chandra 


12989 


2011/03/21 


49.4 


1730 ± 40 


0.31 


Src 


. 8 


XMM-Newton 


0145190201 


2003/02/20 


10.6/24.5/25.2 


7640 ± 90 


2.40 






XMM-Newton 


0145190101 


2003/02/28 


9.4/18.8/19.1 


5460 ± 80 


2.41 


Src 


. 9 


Chandra 


553 


1999/10/19 


5.9 


14 ±4 


4.83 






Chandra 


1454 


1999/10/19 


11.4 


29 ±5 


4.87 






Chandra 


1666 


2001/08/30 


48.6 


120 ± 10 


2.29 






XMM-Newton 


0112980101 


2002/09/28 


13.1/14.7/14.7 


130 ± 20 


4.11 






XMM-Newton 


0112980401 


2002/09/30 


9.0/10.4/10.4 


70 ±10 


4.15 






XMM-Newton 


0112980501 


2002/10/04 


7.5/8.4/8.4 


60 ±10 


4.12 






Chandra 


7698 


2007/06/13 


5.1 


8±3 


4.07 


Src. 


10 


XMM-Newton 


0025541001 


2001/06/19 


6.0/9.1/9.2 


1270 ± 40 


0.55 






Chandra 


11230 


2009/08/11 


24.8 


160 ± 10 


2.30 






Chandra 


10120 


2009/10/24 


10.0 


56 ±8 


1.21 



Notes: a Observation identifiers. b The sum of the good time intervals for each detector. For XMM-Newton observations it is shown in the 
order PN/MOS1/MOS2, a dash indicates that the source was not detected. c Total background-subtracted source counts accumulated in 
the observation (sum of all 3 EPIC detectors for XMM-Newton observations). d Angular separation between the on-axis position of the 
observation and the 2XMM source position. e The Chandra observation of Src. 3 was able to resolve the candidate HLX into multiple 
sources (Figure [TJ. * NGC 4874 is the central galaxy in the Coma cluster, and the Coma cluster has been extensively observed by both 
Chandra and XMM-Newton (in many cases for calibration purposes). There are a number of similar observations from both observatories 
that are not included in this analysis, due to the short exposure time (< 10 ks), but are available in the archives. 



iraf. The specific tools used were: GSFLAT to extract the fiat 
field; GSREDUCE to reduce both the target and calibration 
arc exposures; GSWAVELENGTH to establish the wavelength 
calibration, then GSTRANSFORM to apply it to the exposures; 
GSSKYSUB to subtract the sky background; GSCRREJ to reject 
cosmic ray events; GSEXTRACT to extract one dimensional 
spectra of the source, which were then combined using the 
tool GEMARITH. Gaussian line profiles were fitted where ap- 
propriate to the reduced optical spectrum of the counterpart 
using the CURVEFIT function in IDL version 6.2. 



3 DATA ANALYSIS & RESULTS 

3.1 Spatial analyses 

3.1.1 High resolution X-ray imaging 

As an initial step, we investigated the spatial character- 
istics of the X-ray source detections. One particular issue 
for the XMM-Newton detections of these objects, especially 
given their relatively large distances, is whether we are see- 
ing a single point-like object, or whether multiple objects 
and/or spatially extended components are present within 
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Table 3. Details of archival HST observations, and Galactic ex- 
tinction in the direction of the targets 

Source E(B — V) Prop. ID Instrument Filter 



Src. 3 



Src. 



Src. 8 



0.1 



0.009 



0.011 



6355 WFPC2 F300W 

8597 WFPC2 F606W 

10877 WFPC2 F555W 

F814W 

10861 ACS/WFC F475W 
F814W 



6092 



WFPC2 F450W 
F814W 



Src. 10 0.112 



11575 ACS/WFC F555W 
F814W 



the XMM-Newton beam. Chandra observations offer an ex- 
cellent means of settling this issue, with the much smaller 
Chandra beam able to resolve objects that are confused 
by XMM-Newton, even in regions of high diffuse emission 
(e.g. in the centre of the starb urst galaxies NGC 3256 & 
NGC 3310. IJenkins et alJlJOOJ ). We therefore began by ex- 
amining the deepest Chandra detection for each object in 
our sample, in order to search for resolved structure. 

Only one XMM-Newton detection was resolved to reveal 
finer structure by Chandra, namely Src. 3. This object had 
previously been ident ified as a HLX candi d ate in the spiral 
galaxy NGC 2276 bv lDavis fc Mushotzkvl (J2004T I using the 
XMM-Newton data (Figure [l] top). However, their assumed 
distance for NGC 2276 of 45.7 Mpc is somewhat discrepant 
compared to more recent measurements and overestimates 
the distance to NGC 2276 (cf. Table[T]). This in turn inflates 
the measured lumi nosity. The assum ptions in this paper, and 
the redshift from lde Vaucouleurs et al.l l|1991l ). lead to a dis- 
tance of 33.3 Mpc and hence a lower (although still extreme) 
2XMM luminosity estimate of ~ 6xl0 40 erg s" 1 . In the later 
Chandra observation of this source, the ULX candidate was 
shown to be resolved into three point sources (Figure [TJ bot- 
tom; see also lWolter et alJl2011f ). although only the southern 
two lie coincident with the measured XMM-Newton source 
centroid position. Remarkably, each of the point sources is 
sufficiently luminous to be classified as a ULX, if at the dis- 
tance of NGC 2276, with £ x > 5 x 10 39 erg s -1 . To find 
3 such luminous ULXs within a small region of a galaxy - 
the projected separation of the 3 objects is < 150 pc - ap- 
pears very remarkable, although this should be tempered by 
the realisation that the western regions of NGC 2276 have 
heightened star formation, with a galax y- wide star forma- 
tion rate estimated at 5 — 15 Mq yr^ 1 fWblter et al.ll201ll 
and references therein) . This is plausibly the result of grav- 
itational interaction with NGC 2300 and/or ram pressure 
interaction as NGC 2276 transits through t he bright intra- 
group medium of the NGC 2300 group (e.g. iGruendl et al.l 
11993 ). 

We then used an additional diagnostic test for exten- 
sion on the other, apparently point-like objects. The CIAO 
SRCEXTENT script quantifies the size of each PSF, which 
it does by effectively calculating an average radius for the 
PSF using a Mexican Hat optimisation algorithm. We show 




Figure 1. (top) True colour XMM-Newton image of Src. 3, where 
the rgb colours correspond to the 0.3-1.5, 1.5-2.5 &: 2.5-10 keV 
energy bands respectively. The region circled is the effective 90 
per cent uncertainty in the XMM-Newton position of the source. 
(bottom) Equivalent Chandra image. The XMM-Newton detec- 
tion is resolved into three separate objects, two of which border 
on the 90 per cent XMM-Newton error region. They are referred 
to as labelled here throughout the rest of the paper. These and 
all subsequent images are aligned such that north is up. 



the results of this in Table [4] These results are consistent 
with expectations for the Chandra PSF for both on-axis and 
slightly off-axis objects, where an extent of ~ 0.5 arcseconds 
is equivalent to the 75 per cent encircled energy radius for 
1.5 keV photons (see Fig. 4.6 of the Chandra Proposer's 
Observatory Guide; also the Chandra High Resolution Mir- 
ror Array calibration documentation 11 !). We have verified 
the off-axis extents are consistent with point-like objects 
by comparison to the population o f point-l ike ULXs in the 
galaxy NGC 4490 (seelRoberts et al.ll2002l;[Fridriksson et all 
120081 ; (Gladstone fc Robertdl2009l : lYoshida et al.ll2010r i. eight 
objects that lie up to 3 arcminutes off-axis in three separate 
Chandra ACIS-S observations. Our extents are consistent 
with the SRCEXTENT-measured extents for these ULXs over 
the full range of off-axis angles seen. The only outlier is Src. 
9, but even this is consistent at the 90 per cent confidence 
level. We conclude that there is no evidence for significant 
extended emission for any object, meaning that each ULX 
candidate remains consistent with a single, point-like object 
at the high X-ray spatial resolution of Chandra. 



11 http://cxc.harvard.edu/proposer/POG , 
http : //exc . harvard . edu/cal/Hrma/Index . html 
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Figure 2. DSS2 images of the habitats of the 10 objects in the selected sample. The ordering of panels matches the numbering of the 
sources, starting at top left with Src. 1 and reading across each row. Data from the DSS2-blue survey is shown where available, but in the 
two cases this was not possible (Src. 2 & Src. 7). The DSS2-red filter data is shown instead. Each image is 5 arcminutes across, aligned 
such that North is up, and centred on the nucleus of the host galaxy. The greyscale is arbitrary in each panel, and displayed such that 
the internal structure of each galaxy can be seen. The position of each object is marked by a circle, except where HST data is available 
where instead a box is shown matching the size of the image in Fig. [3] Note that the actual position errors are considerably smaller than 
the size of the circle shown. 
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Table 4. Chandra source characteristics 



Table 5. HST photometry for candidate counterparts 



Source 


Obs ID 


RA, dec 


Q a 


Extent 








(J2000) 


(arcmin) 


(arcsec) 


Src. 


1 


10562 


011942.8 + 032422 


0.29 


0.48 ± 0.03 


Src. 


2 


12988 


024025.6 - 082430 


0.32 


49+0-02 
u -^ a -o.oi 


Src. 


3a 


4968 


072648.1 + 854554 


1.13 


0.49 ± 0.04 


Src. 


3b 


4968 


072647.9 + 854552 


1.10 


_b 


Src. 


3c 


4968 


072648.3 + 854549 


1.06 


45+0- 07 


Src. 


4 


12990 


120405.8 + 201345 


0.30 


0.50 + 0.1 


Src. 


5 


7863 


121856.1 + 142419 


1.84 


0.55 + 0.05 


Src. 


6 


10672 


125939.8 + 275716 


1.93 


62+ 01 
u ' DZ -0.09 


Src. 


7 


12989 


134404.2 - 271411 


0.31 


0.50 + 0.02 


Src. 


8 C 


- 


- 


- 


- 


Src. 


9 


1666 


163614.1 + 661410 


2.29 


0.7 + 0.1 


Src. 


10 


11230 


230457.7 + 122029 


2.30 


0.54 + 0.07 



Source ID Filter 



'/ 



'f 



Notes: Source characteristics are measured from the Chandra ob- 
servation with the greatest number of counts accumulated from 
the ULX candidate in question. a Off-axis angle of source. b 
SRCEXTENT was unable to constrain the size of this point-like 
object due to the proximity of the brighter Src. 3a. However, 
inspection shows its PSF appears very similar to its two verified 
point-like neighbours. c Src. 8 has not previously been observed 
with Chandra. However, a Chandra observation of the source is 
scheduled. 

3.1.2 Source locations and optical counterparts 

A second benefit of obtaining Chandra data for the objects 
in our sample is the accuracy of its astrometry, with sub- 
arcsecond positioning expected even without applying as- 
trometric corrections. We were therefore able to derive ac- 
curate positions (shown in Table [4} and so examine multi- 
wavelength data in order to search for counterparts, and 
identify the environments hosting these objects. Initially we 
investigated the locations of the objects in our sample by 
overlaying their positions on DSS2 images, which we show 
here as Fig. [2] Seven objects are associated with spiral galax- 
ies; in the six cases where the galaxy is close to face-on, the 
ULX candidate appears to be located on or close to the spi- 
ral arms, in most cases in the outer regions of the disc. In the 
one edge-on case (NGC 5907, the galaxy containing Src. 8) 
the source is located in the plane of the disc, with a projected 
location coincident with a strong dust lane. We also searched 
for catalogued objects near to the positions of our objects 
using the NED and SIMBAD databases^- In the case of 
the spiral-hosted objects, three were in close proximity to 
catalogued HII regions: Src. 1 i s ~ 2 arcseconds from HII 
region 46 of lHodee fc KennicuttJ (119831): Src. 3 is within ~ 6 
arcseconds of regions 61 & 65 of lHodee fc KennicuttJ |l983j); 
and Src. 10 i s with in ~ 2 arcseconds of regions 197 & 898 of 
iRozas et al.l |l999j). These associations imply a strong pos- 
sibility that these objects are associated with ongoing star 
form ation, as appears the case for m any less luminous ULXs 
(e.g. ISwartz. Tennant fc Sori3l2009h . 

The three objects associated with elliptical galaxies 
present a more diverse subsample. As Fig. [2] shows, the ob- 
ject near to NGC 4065 (Src. 4) appears to have a faint, 
point-like counterpart. Subsequent to the selection of the 

12 http : //ned . ipac . caltech . edu , http : //simbad . u-strasbg . f r 



Src. 3a 



Src. 3b 



Src. 3c 



Src. 6 



Src. 



Src. 10 



F300W 22 + 4 -11.1 

F555W 23.1 + 0.9 -9.8 

F606W 22.9 + 0.8 -10.0 

F814W 21.8 + 0.8 -11.0 



F300W 


- 


- 


F555W 


25 + 2 


-7.9 


F606W 


25 + 3 


-7.9 


F814W 


23 + 5 


-9.8 


F300W 


_ 


_ 


F555W 


- 


- 


F606W 


24 + 3 


-8.9 


F814W 


24 + 5 


-8.8 


F475W 


22.0 + 0.2 


-13.0 


F814W 


21.2 + 0.2 


-13.8 


F450W 


21.5 + 0.4 


-9.4 


F814W 


- 


- 


F555W 


26 + 2 


-6.9 


F814W 


26 + 5 


-6.8 



Notes: a Apparent magnitudes, for the filter/, in the Vega photo- 
metric system. b Absolute magnitudes, corrected for foreground 
extinction as per Table [3] 

sample this object was revealed to be a background QSO 
by SDSS spectroscopy (object SDSS J120405. 83+201345.0). 
We therefore do not consider it further in the main body of 
this paper, although we briefly report its properties in Ap- 
pendix A. The second object associated with an elliptical, 
Src. 6, is actually found to be in the close proximity of a 
satellite galaxy (SDSS J125939.65+275714.0) of the central 
dominant galaxy of the Coma cluster, although it is dis- 
placed from the satellite galaxy nucleus by ~ 4 arcseconds. 
The final object, Src. 7, is the most X-ray luminous object 
in the sample, although also amongst the most distant. Its 
host is a pe culiar SO object that appears to be crossed by 
a dust lane (|Laubertsl 1 1982T ) . However, there is no obvious 
counterpart to this object in the DSS image. 

A minority of the objects in the sample are also cov- 
ered by HST imaging observations, as listed in Table [3] 
We show these four datasets in Fig [3] where we display 
the immediate region covering the candidate ULX position 
at the full HST spatial resolution. We overlay an X-ray 
position error circle on top of the image to indicate the 
likely position of each ULX candidate. This is calculated 
as the combination of position errors from the astromet- 
ric uncertainty of the Chandra data, the X-ray source cen- 
troi<L_and the astrometric registration of the HST data, as 



pe nRoberts. Levan fc Goadl (|2008T ). Possible point-like can- 
didate counterparts are apparent within the X-ray position 
error circle for all of the objects. Aperture photometry was 
performed on the brightest identified counterpart in each 
case (which in most cases is the only counterpart, although 
there are possibly fainter objects within the crowded regions 
in NGC 2276), as described in Section 2.4, and we list the 
derived magnitudes of the counterparts in Table [5] 

The brightest individual objects within the error cir- 
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Figure 3. HST images of the regions containing four of the candidate ULXs. Top left — Src. 3a, b & c. This HST data is displayed as 
a 3-colour image, with the F814W filter providing the red data, F555W the blue, and green taken as an average of the two. The data is 
displayed on a linear scale to aid the visualisation of the counterparts. Top right - Src. 6. The image is 3-colour, with the F814W data 
in red, F475W in blue, and an average of the two filters in green. The image is displayed logarithmically to permit the separation of the 
counterpart from the galaxy to be seen. Bottom left - Src. 8. A 3-colour image using F814W as red, F450W as blue, and an average of 
the two as green. This is scaled linearly for display purposes. Bottom right - Src. 10. A single, linearly greyscaled image in the F555W 
band is presented to emphasize the faint counterpart. 



cles for the ULXs in NGC 2276 (Src. 3a-c) display ab- 
solute magnitudes in the ranges Mf606\v ~ —8 — ► —10 
and Mf8i4\v ~ —9 — ¥ — 11, albeit with large uncertainties. 
Archival HST F300W and F555W observations covering the 
field of Src. 3a-c are also available, but the magnitude of 
the brightest optical source within the error circle of Src. 
3c is not constrained by the F555W observation, and simi- 
larly for both Src. 3b and c in the F300W observation. Such 
magnitudes are more luminous than typical individual su- 
pergiant stars, although examples of luminous blue variable 
stars with similar magnitudes are kno wn (e.g. Cyg OB2-304 
with Mv ~ — 10.6; [Massev et alj200lT ). However, these mag- 
nitudes are also reasonable for young massive star clusters, 
such as one might find in an intense star forming region, and 
even somewhat unremarkable - for example more than 50 
young clusters in the Antennae have ab solute magnitudes 
of M v < -11.3 IWhitmore et al.lll999T ). Src. 6 is a very 
interesting case, as it is located in an early type galaxy. 
However, the possibility of a globular cluster host seems 
rather unlikely, as its absolute magnitude is larger than 
would be expected for such an object. A very interesting 



alternative is that it could be located in an Ultra-Compact 
Dwarf (UCD) galaxy, a number of which are also known to 
be located within t he central regions of the Coma cluster 
i|Madrid et al.ll2010T l. Given that UCDs are thought to fill 
an evolutionary step between the largest globular clusters 
and dwarf spheroidal galaxies, such an environment might 
be an ideal place to look for IMBHs. However, the colour 
of the object is anomalously blue (F475W-F814W w 0.8) 
compared to other UCDs in the same field, which typically 
have F475W— F814W colours of ~ 1.5 to 2, and so is per- 
haps more indicative of a background QSO. The next ob- 
ject, Src. 8, is a very blue object, not being detected in the 
F814W filter, and hence must be located at the near side 
of the dust lane in its host galaxy (NGC 5907) . Its absolute 
magnitude of ~ —9.5 in the F450W filter is too large for a 
single stellar object, but is potentially indicativ e of a young, 
blue stellar cluster fe.g. IWhitmore et al.lll999f ). Finally, the 
counterpart to Src. 10 appears consistent with a single giant 
star, although again its magnitudes are poorly constrained 
by the data and so no further constraints from its colour 
information are forthcoming. 
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3.1.3 Optical spectroscopy of a counterpart 

Given the interesting nature of the counterparts to these 
objects, further follow-up observations to identify their na- 
ture are critical. To this end, we have already obtained a 
follow-up low-dispersion optical spectrum of the possible 
UCD galaxy counterpart to Src. 6 using the GMOS long 
slit spectrometer on the Gemini-North telescope. Unfortu- 
nately, in this case the spectrum revealed emission lines con- 
sistent with a relatively high redshift QSO (z m 3.25), thus 
explaining its blue colours. We therefore identify this object 
as the second contaminant in our sample, and do not con- 
sider it further in the analysis (although, we again describe 
its properties in Appendix A) . 



3.1.4 Background contamination estimate 

Out of the ten objects in our initial sample we have there- 
fore identified two background contaminants. But how many 
should we expect? To answer this question we have revisited 
the background contamination calculations presented for the 
parent 2XMM-DR1/RC3 sample in IWalton et al.l (2011b). 
We repeat the process presented in that paper, with appro- 
priate changes for consistency with this sample, i.e., we set 
the luminosity limit to be 5 x 10 40 erg s _1 , and only con- 
sider galaxies within 100 Mpc. From this we predict a false 
detection rate of ~ three objects with apparent luminosi- 
ties above 5 x 10 40 erg s _1 across our entire sample, albeit 
with large uncertai nties stemmin g from the small number 
statistics. As with IWalton et al.l (2011b), the background 
contamination is expected to be higher for elliptical galax- 
ies, where we predict ~ two QSO contaminants, compared 
to ~ one for spiral galaxies. This compares favourably to 
what we have detected; the two identified QSOs are indeed 
located close to elliptical galaxies, matching the prediction. 
On the other hand, the low predicted background rate for 
spiral galaxies (one out of seven detections) provides fur- 
ther confidence that the majority of these objects are real 
luminous ULXs, associated with their parent galaxies. 



3.2 X-ray spectral analysis 

Where sufficient counts were available to enable X-ray spec- 
tral analysis, the spectra were fitted with basic spectral mod- 
els using XSPEC v. 12.6.0. Observations with > 250 counts 
were grouped using the FTOOL GRPPHA to a minimum of 
twenty counts per energy bin, to produce bins with Gaus- 
sian errors and hence permit the use of \ 2 statistics in fitting 
the data. Where only 100 — 250 counts were available, models 
were fitted t o ungrouped X-ray spectra using Cash statistics 
i|Cashlll979T ). The data were not background subtracted for 
the Cash analyses; instead in each of these cases the back- 
ground spectrum was first constrained from the background 
aperture, and then added to the fit as a constant component, 
normalised for the relative sizes of the extraction regions. In 
all cases the fitting was limited to the 0.3-10 keV band. 

The spectra were fitted with both an absorbed power- 
law continuum and an absorbed multi-colour disc blackbody 
(diskbb in XSPEC, hereafter MCD) models, with the latter 
representative of the emission from a stan dard optically- 
thick, geometrically thin accretion disc (|Mitsuda et al.l 



1 1984 ) . Absorption w as modelled by the tbab s mod el, us- 
ing abundances from IWilms. Allen fc McCravl (J2000). Two 
absorption components were included, the first of which 
was fixed to the Galactic line-of-sight column to the object 
i|Dickev fe Lockmarjl990l . cf. Table [TJ, and the second left as 
a free parameter to model intrinsic absorption in the source 
and/or its host galaxy. Where multiple XMM-Newton EPIC 
detectors were operating, an additional multiplicative con- 
stant was included in the spectral model. This was fixed 
to 1.0 for the PN detector (or MOS1 where no PN data 
were available) and varied freely for the others, in order to 
account for the calibration uncertainties between the EPIC 
detectors. In all cases with equally sized extraction apertures 
the values of this parameter remained within a few percent 
of unity, and the values scaled appropriately where smaller 
extraction regions were necessary on individual detectors. 

Table [S] shows the parameters resulting from fitting the 
two simple models to the source data. Plots of the absorbed 
power-law spectral fits, unfolded from the detector response, 
are displayed in Fig. [3] for all observations with greater than 
250 source counts. An absorbed power-law continuum spec- 
trum provides a statistically acceptable fit (i.e., null hypoth- 
esis probability > 5 per cent, so rejection likelihood < 95 
per cent) to data from all but two observations. These ex- 
ceptions are the least luminous of the resolved source triplet, 
Src. 3c, which was not well characterised by either model; 
and a marginal rejection of the latter, much less luminous 
observation of Src. 10, which was the sole dataset acceptably 
fitted by the MCD model alone. In contrast, eight spectra 
were acceptably fitted by the absorbed power-law model, but 
rejected a MCD fit at greater than the 95 per cent level, and 
11 spectra were acceptably fitted by both models. In all, the 
absorbed power-law fit was statistically preferred (i.e., pro- 
vided the closer fit to the data) in 15/21 cases, compared to 
5/21 for the MCD, although many of these preferences were 
marginal in nature. 

The large majority of the spectra that were acceptably 
characterised as absorbed power-law continua all had rather 
hard photon indexes, with F = 1.4-2.2. The absorption col- 
umn densities for these fits were mostly moderate (where 
well constrained), with values in the range ~ 1 — 3 x 10 21 
cm -2 , although the column for Src. 8 was notably higher 
at ~ 10 22 cm" 2 , which may be due to its location in an 
edge-on system. It should be noted however, that such an 
interpretation may bring into doubt any physical associa- 
tion with the particularly blue potential counterpart to this 
source. We discuss this further in section 4.1. The spectra 
that were acceptably fitted by the absorbed MCD model 
had - with one notable exception - inner disc temperatures 
of fcTi n ~ 1.0-1.9 keV, and rather low modelled absorp- 
tion columns at < 1 x 10 21 cm -2 , where well-constrained. 
The exception was the fainter observation of Src. 10, where 
the MCD model was the only acceptable fit, and the spec- 
trum was remarkably soft with an inner disc temperature 
fcT in ~ 0.4 keV. 

An inner disc temperature in the range 1-2 keV 
is fairly typical for a sub-Eddington stellar mass black 
hole in the thermal dominant state, with the accre- 
tion disc extending to the innerm ost stable circular or- 
bit l|McClintock fc Remillardl l2006h . However, for objects 
at close to Lx ~ 10 41 erg s _1 it would appear some- 
what anomalous, and certainly not indicative of the pres- 
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Figure 4. XMM-Newton EPIC PN (black), MOS1 (red), MOS2 (green) or Chandra spectra for observations with greater than 250 source 
counts, all of which are unfolded from the detector response. Multiple observations of sources are shown side-by-side where possible. 
Where there are two or more observations of a source with the same detector, plots are shown on the same scale for direct comparison. 
Over-plotted lines show the best fitting absorbed power-law spectra (Table |6|, with one absorption component set to the Galactic value 
and another free to model intrinsic absorption in the ULX candidates. A range of photon indexes between 1.4-2.2 were obtained. 
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Figure [4) (continued) 
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Table 6. X-ray spectral modelling 









Absorbed power- 


-law 


Absorbed MCD 




Source ID 


Obs ID 


N H a 


V b 


X 2 /dof c 


N H a 


kT in d 


X 2 /dof c 


/x,PL 


Src. 


1 


0200780101 


o.26±°:i 8 


2.2 ±0.2 


54.8/43 


< 0.09 


1.0 ±0.1 


34.5/43 


o n+°- 4 

S - U -0.8 






10562 


14+ - 07 


1.8 ±0.2 


27.9/25 


< 0.03 


i O+0.2 
^-0.1 


[39.9/25] 


1 99+0- 08 






0601670101 


0.10 ±0.02 


1.80 ±0.08 


132.3/131 


< 0.004 


1 22+ 01 
± - zz -0.08 


[227.8/131] 


1 4Q+0-05 
i -* M -0.05 


Src. 


2 


0093630101 


20+ 006 
u - zu -0.05 


2.2 ±0.2 


43.3/41 


< 0.04 


0.9 ±0.1 


[74.4/41] 


6.6±g;J 






0306230101 


n -,0+0.04 
U - 18 -0.03 


1.61 ±0.09 


126.6/133 


0.03 ± 0.02 


1.6 ±0.1 


154.7/133 


s.4±S:§ 






0553300401 


0.16 ±0.03 


1 kq+0.08 
1 - oy -0.07 


219.6/214 


02+ 002 
uuz -0.01 


1.6 ±0.1 


[281.6/214] 


10.8±g-2 






12988 


0.26 ±0.04 


1 Q+01 

L - i '-o.o9 


122.4/132 


0.04 ± 0.03 


1 o + O.l 

^-o.og 


152.6/132 


s.o±8:S 


Src. 


3 


0022340201 


0.10 ±0.04 


-, 42 +0.09 


107.1/92 


< 0.01 


I Q+02 

J^-o.i 


104.8/92 


4 6+ ' 2 
4 - D -0.2 


Src. 


3a 


4968 


5 +0 - 4 
U.O_ 2 


1.9 ±0.3 


16.4/12 


< 0.4 


i 9+O.5 
L - z -0.2 


15.1/12 


0.76+°'« 6 


Src. 


3b 


4968 


0.4 ±0.2 


9 1+0.3 
z ' i -0.4 


10.9%' 


< 0.2 


i O+0.4 
1 -°-0.2 


13.9%' 


0.35+Of 


Src. 


3c 


4968 


1 5+ 1 


1.7 ±0.6 


[100%'] 


o.9i° :l 


9 "}+ 2 
z - J -0.8 


[100%'] 


0.39+ , , 6 


Src. 


5 


0147610101 


26+ ' 06 
u ' ZD -0.05 


1.5 ±0.1 


74.1/86 


omtVol 


1.8 ±0.2 


94.1/86 


4.4±g;l 






7863 


< 0.5 


1 4+0-5 


6.8/10 


< 0.2 


i 7 +0.9 
i -'-0.4 


8.7/10 


6.8+0- 8 


Src. 


7 


0092970801 


< 0.2 


1.7 ±0.3 


15.0/16 


< 0.05 


1.2 ±0.2 


16.7/16 


2.8±gi 






12989 


0-0810.05 


1.63t°f 


78.8/69 


< 0.007 


1 34+ 01 


[103.0/69] 


q 9+0.2 
°- z -0.3 


Src. 


8 


0145190201 


0-97±° ° l 


1 SS+0' 07 


409.5/366 


0.54 ± 0.04 


1 66+ - 08 


[420.0/366] 


15.8185 






0145190101 


n O-+0.07 


1.68 ±0.07 


267.0/251 


0.49 ± 0.04 


1.9 ±0.1 


258.5/251 


14.718:5 


Src. 


9 


1666 


< 0.2 


1 6 +0 ' 4 


1.1%' 


< 0.06 


1 4+0-5 


[98.4%'] 


n S4+ - 06 

u - J4 -0.09 






0112980101 


< 0.2 


1 8+0- 4 


83.3%-' 


< 0.03 


1 5+°- 4 


[97.8%'] 


nfil +0.06 
u - oi -0.09 


Src. 


10 


0025541001 


n 97 +0.08 
u - z '-0.07 


1 Q+0-2 

J-^-o.i 


61.3/62 


< 0.07 


1 4+0-2 


66.1/62 


4.7lg;I 






11230 


0.6 ±0.2 


4 2+ - 9 
4 - Z -0.8 


[96.7%'] 


< 0.3 


40+ 01 


57.3%' 


o.28t8:° 2 4 



Notes: Errors shown for spectral parameters are 90 per cent confidence ranges. a Absorption column density external to our Galaxy 
(x 10 22 cm -2 ). b Power-law photon index. c Statistical goodness of the fit, in terms of the x 2 statistic and number of degrees of freedom. 
Numbers in bold indicate which model provided the best fit to the data, whilst numbers in square brackets indicate a fit that is rejected 
at 2cr significance. d Inner-disc temperature (keV). e Observed 0.3-10 keV flux, in units of 10 -13 erg cm~ 2 s -1 , base d on the power-law 
fit. * Total counts from this observation were insufficient for x 2 fitting, instead Cash statistics were used IICasnlll979r ). Goodness of fit is 
shown instead of x 2 /dof. 



ence of IMBHs, which should have much cooler disc- 
dominated spectra (e. g. iKaaret et al.l 12003 : iMiller et al.l 
120031 : iMiller et al.ll2004l ). as indeed does seem to be the case 
for the fainter dataset for Src. 10. We investigated whether 
there was any evidence for cool disc components in the other 
spectra by fitting a two-component, absorbed MCD plus 
power-law model to the best datasets (the 12 spectra with 
> 1000 counts), where we had the best opportunity of de- 
convolving the spectra into multiple components. In all but 
two cases there was no discernable improvement to the ab- 
sorbed power-law spectral fit with the addition of a cool disc 
component. Two cases offered marginal improvements: the 
fit to the original XMM-Newton detection of Src. 3 (before 
being resolved by Chandra) improved by Ax 2 = 10.9 for 
two additional degrees of freedom with a ~ 0.1 keV disc; 
and the first, most luminous detection of Src. 1 improved 
by Ax 2 = 13.1 for two additional degrees of freedom with 
a ~ 0.07 keV disc. However, in the latter case the MCD 
model alone still offered the most acceptable fit to the data. 
This lack of strong cool disc components in the sample is 
perhaps unsurprising given the lack of strong soft residuals 
in the spectra presented in Fig. [4] 

However, one object stands out as possessing interest- 
ing residuals in Fig. [4] Both observations of Src. 8, the 
highest quality data in the sample with > 5000 counts per 
dataset, appear to show a high energy spectral turnover be- 
yond ~ 5 keV. Interestingly, this feature is the key diagnos- 
tic of the putative ultraluminous state, that is seen in many 



less-luminous ULX spectra (jGladstone et al.ll2009l ). We note 
that the other main diagnostic of this putative state, a soft 
excess, is not seen; although this may be due to the object 
being buried in an edge-on galaxy, and hence its soft X- 
rays being relatively highly obscured, as is suggested by its 
spectral fits (cf. Table [5]). Interestingl y, this is very compa- 
rable to the object NGC 4517 ULX1 J Walton et all 2011a), 
which is also located in an edge-on galaxy and has a highly 
absorbed soft spectrum, with a break at high energy. Confir- 
mation was sought o f the significanc e of th e spectral break 
using the method of IStobbart et al.l (|2006t ). i.e., by fitting 
unabsorbed power-law and broken power-law models to the 
2-10 keV X-ray spectra of both observations of Src. 8. A 
broken power-law model provided a better fit to both obser- 
vations, with very significant F-test probabilities of 1 x 10 -12 
(0145190201) and 4 x 10" 8 (0145190101), hence the spec- 
trum is confirmed to be consistent with the ultraluminous 
state. 

Both Fig. [4] and Table [6] demonstrate that the indi- 
vidual objects also displayed both flux and spectral shape 
variations between different observing epochs. In some cases 
the spectral variability was quite subtle, for example a hard- 
ening of the absorbed power-law spectrum with flux in Src. 
2. In other objects the contrasts were much starker, partic- 
ularly in Src. 1 and Src. 10. In the former case the object 
transits from an absorbed disc-like spectrum at high flux 
to an absorbed power-law spectrum in the two subsequent 
spectra where the flux is a factor > 4 lower, and in the latter 
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Table 7. Luminosity and variability of the extreme luminosity 
ULXs 



Obs ID 


Lx a 


F var (200 s) b 


F var (2 ks) b 


Src. 1 








0200780101 


10.2+ / 5 


< 0.2 C 


- 


10562 


2 5+ - 1 


- 


0.11 ±0.06 


0601670101 


1 Q0+ 007 


- 


< 0.1 c 


Src. 2 








0093630101 


2.8 ±0.2 


- 


- 


0306230101 


og+0.1 
O.D_g 2 


- 


0.05 ±0.04 


0553300401 


4 6+ 01 
*'°-Q.2 


0.09 ± 0.04 


0.10 ±0.02 


12988 


3 4+0.1 
°- -0.2 


< 0.1 c 


0.13 ±0.02 


Src. 3 








0022340201 


6.1 ±0.3 


- 


0.08 ±0.04 


4968 - a 
4968 - b 
4968 - c 


1 01+ 008 

J - ux -0.2 

0.46+H 6 
0.521H 8 


- 


- 


Src. 5 








0147610101 


c u+0.2 
°* 8 -0.3 


0.13 ±0.06 


- 


7863 


9 +1 

9 -3 


< 0.3 C 


- 


Src. 7 








0092970801 


30± 2 


< 0.3 C 


- 


12989 


35± 2 


- 


< 0.1 c 


Src. 8 








0145190201 


4.19±°;° 9 


0.07 ±0.02 


- 


0145190101 


3.9 ±0.1 


0.11 ±0.02 


- 


Src. 9 








553 


2.9 ±0.9 


- 


- 


1454 


3.2 ±0.6 


- 


- 


1666 


3.8+r 


- 


- 


0112980101 


e.str 


- 


- 


0112980401 


5.9 ±0.8 


- 


- 


0112980501 


5±1 


- 


- 


7698 


1.9 ±0.7 


- 


- 


Src. 10 








0025541001 


6.1 ±0.3 


0.10 ±0.06 


- 


11230 


0.36+. ,; , 5 


- 


- 


10120 


0.40 ± 0.05 


- 


- 



Notes: a 0.3-10 keV observed luminosity, in units of 10 erg s _1 , 
calculated as per the text. The quoted errors are the la confidence 
ranges on the luminosity. ^Fractional variability amplitude. Val- 
ues and errors are shown where variability is detected at > 1 
a significance (the quoted errors are the la confidence ranges). 
Upper limits shown are 3<r limits, and are calculated where the 
variability was less than white noise. The FTOOL LCSTATS was used 
to provide this upper limit. 



case the spectrum changes from hard and power-law-like at 
high flux to much softer and disc-like after a drop in flux 
of nearly a factor 20. We discuss these changes further in 
Section [5] 



3.3 X-ray timing 

3.3.1 Short term variability 

In addition to extracting X-ray spectra, we also extracted 
background-subtracted light curves for each observation. A 
representative selection of these light curves is shown in 
Fig- El covering the differing data qualities, exposure lengths 
and amounts of variability in the data (see below). The light 



curves are binned such that there are ~ 25 counts per ele- 
ment of temporal resolution. 

The fractional variability amplitude, F val , was used to 
determine whether statistically significant short term vari- 
ability was present in the incident count rate of the sources. 
Here we define F var as 



*■ var — 



M Lisll 1 ' X ) N^i = l a ' 



2 

crr.i 



(1) 



as pe nVaughan et al.l ([20031 ). where N is the number of tem- 
poral bins used, Xi are the individual count rates per bin, 
x is the mean count rate and o~\ rr i are the variances on 
the individual count rates. Variability was tested for on two 
separate fixed timescales, 200 s and 2 ks, and the results 
are shown in Table [7] Unfortunately, in most of the 24 ob- 
servations the low count rates resulted in insufficient counts 
per temporal bin (< 20), or short exposure times limited 
the available number of temporal bins (again, < 20), which 
made such an analysis impractical. These datasets are shown 
as the blank spaces in Table [7] In fact, the data was only 
adequate for these analyses in nine observations at 200 s, 
and seven at 2 ks. In determining _F var temporally incom- 
plete light curve bins (e.g. due to overlapping with the end 
points of GTIs, or the on/off times of instruments) were ex- 
cluded so that the errors on the data points were normally 
distributed. Short term variability was detected at greater 
than 2<t significance in observations of four (and at greater 
than 3<t in two) of the eight sources in the sample, most no- 
tably from both observations of Src. 8 on 200 s timescales. 
The other marginally significant detections and upper limits 
were of sufficient quality to constrain values that were con- 
sistent with similar levels of variability, mainly at ~ 10 per 
cent of the total count rate. 



3.3.2 Long term variability 

In order to study the long term variability of the sources, 
fluxes were obtained for each observation, that were con- 
verted to luminosities using the distances in Table [l] Where 
an observation contained sufficient counts to allow spectral 
fitting, fluxes were determined using the absorbed power- 
law continuum spectral fit. For observations with insufficient 
counts for spectral fitting, light curves were \ 2 fitted with a 
constant count rate, and the resultant best fitting value was 
converted to a flux using PlMMa 13 !. A characteristic spectral 
index of 2.0 and an intrinsic absorption column of 1.0 x 10 21 
cm -2 were assumed, except for observation 10120 of Src. 10, 
where the parameters of the much softer disc-like spectrum 
seen in observation 11230 were used instead. The observed 
luminosities are shown in Table alongside the F vaT values. 
Most objects display long term luminosity variations at 
some level between the observations. In three objects the ini- 
tial high luminosity that led to them being included in the 
sample is not maintained in subsequent observations, and 
appears to drop by factors > 4. Indeed, as a result of spec- 
tral fitting, two objects appear to be slightly under-luminous 
compared to the selection criteria of the sample (Src. 2 & 
Src. 8), although in both cases their peak luminosity is still 
considerable at ~ 4 x 10 40 erg s _1 . In contrast, the most 
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Figure 5. A selection of background-subtracted light curves of the ULX candidates, chosen to show a variety of behaviours and data 
quality: continuous, moderate count rate data (top left & top right); sparse data (middle left & middle right); low count rate data (bottom 
left); variability about a mean count rate (top left & top right) and an apparent maintained shift in the mean count rate (bottom left 
& bottom right). XMM-Newton light curves are summed over all EPIC detectors, except for 0093630101, where only PN and MOS2 
detections were used. The dashed line shows the minimised \ 2 best fitting constant count rate for each observation. 



luminous object, Src. 7, remains persistently extremely lu- 
minous at ~ 3 x 10 41 erg s _1 in both observations obtained 
to date. 



4 COMPARING THE CHARACTERISTICS OF 
THE SAMPLE SOURCES TO OTHER ULXS 

In the previous sections we have examined the characteris- 
tics of a small sample of the most luminous ULXs. Here, we 
will attempt to place these results in the wider context of 
ULX and HLX studies. 



4.1 Spatial properties 

Similarly to other ULXs, all but one of the sample re- 
main unresolved and point-like at the X-ray spatial reso- 
lution of the Chandra observatory, which is the highest cur- 
rently available. In the one case where the original XMM- 



Newton source detection was resolved to reveal underlying 
structure, Src. 3, it resolved into three luminous ULXs, each 
of which is point-like. However, the sum of luminosities in 
the three resolved objects is substantially below the lumi- 
nosity of the original XMM-Newton detection. Moreover, by 
virtue of the short-term variability in the original detection 
we can conclude that it is likely that just one of the three ob- 
jects was substantially brighter in the first epoch, such that 
it was the dominant source of flux in the XMM-Newton de- 
tection. As it is unclear which of the three subsequently 
resolved objects was brightest in the original data, we re- 
tain all three in this work. However, for the purposes of 
constraining the relative changes within this object (e.g. its 
luminosity with time), we base the second epoch measure- 
ment on the brightest of the three resolved objects, Src. 3a, 
in the remainder of this paper, which effectively provides a 
lower limit on the long term change in flux. 

The locations of the candidate ULXs within their host 
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galaxies may also be used as a comparator. The objects 
that lie coincident with face-on spiral galaxies all appear 
associated with the outer regions of their hosts, close to 
spiral arm structure. In three galaxies that have had their 
Hn regions mapped, the ULXs all lie close to known Hn 
regions. This evidence would seem indicative of these ob- 
jects being associated with recent star formation; and this 
appears in common with less luminous ULXs. Indeed, an 
association with star formation has long been qualita- 
tively established by the discovery of populations of ULXs 
in star forming galaxies in the early years of the Chan- 
dra mission (e.g . iFabbiano et al.l l200ll : iRoberts et al.l 120021 : 
iGao et al.1 120031 ) . with subsequent studies quantifying the 
relat ionship with the amount of ongoing star formation 
(e.g. ISwartz et alJlJOOJ : iGilfanov. Grimm fc Sunvaevll20o4 
ISwartz et al.l 12003 ). Furthermore, the quoted separation of 
the extreme ULXs from known Hn regions is of the or- 
der ~ 200 pc in NGC 470 and NGC 7479, consistent with 
the separation of ULX s from stellar clusters in the Anten- 
nae (IZezas et al. 20021) and other nearby starforming galax- 
ies (jKaaret et al.l 12004! ) . In the two face-on galaxies with 
HST imaging, NGC 2276 and NGC 7479 (both ~ 33 Mpc 
distant, where the ~ 1 arcsecond diameter error region 
equates to ~ 160 pc), the objects appear near to likely star 
forming complexes in the spiral arms. All four of these ULXs 
possess candidate counterparts, consistent with stellar clus- 
ters in the three NGC 2276 ULXs (Src. 3a-c), and a possi- 
ble supergiant stellar counterpart to the NGC 7479 object 
(Src. 10). Thus the objects remain at the very least associ- 
ated with a star forming environment, even in HST imag- 
ing. This is also a characteristic in common with two of the 
three best studied HLXs; M82 X-l is lo cated in the cen- 
tre o f an archety pal starburst galaxy (e.g. lMatsumoto et al.l 
l200ll ; iKaaret et al.l [200lh . and object N10 is in the famou s 
starforming ring of the Cartwheel galaxy (|Gao et al.M 2003). 
Src. 8 presents an interesting case; it lies coincident with 
a dust lane in the edge-on galaxy NGC 5907, and is rela- 
tively highly absorbed in the X-rays, consistent with being 
embedded within the dust lane; yet has a very blue counter- 
part, that is also relatively luminous with Mfasow = —9.4. 
It is difficult to interpret this as a young stellar cluster or 
a single giant star as there is no evidence for further star 
formation in its vicinity; but this object cannot be a back- 
ground object as the high extinction through the NGC 5907 
dust lanes would not permit a background object with blue 
colours to be visible. Given the X-ray characteristics of this 
object discussed later, it would seem to be a plausible ULX; 
therefore its blue colour and high absolute magnitude sug- 
gest its optical emission is dominated by reprocessed X-ray 
emission from the outer regions of the accretion disc, as 
suggested for some fainter ULX o ptical counterparts (e.g. 
iTao et al.ll201l| ; iGrise et alj l2012h and ESO 243-49 HLX- 
1 (|Farrellet alll2012J ; ISoria et al.ll2012i ). although the high 



X-ray absorption column remains a problem for this inter- 
pretation. 

After the exclusion of two of the sample near elliptical 
galaxies as background QSOs, we are left with Src. 7 alone in 
this type of habitat. It is somewhat surprising that we find 
any HLX in this type of environment, given their lack of star 
formation. Indeed, previous results show that the X-ray lu- 
minosity functi ons of elliptica l galaxies ar e steep (e.g. com- 
plete sample of IWalton et al.l 2011b; also iKim fc Fabbiand 



I2010T) . and llrwin. Bregman fc Athevl (J2004T ) demonstrated 
that the numbers of ULX candidates with luminosities above 
2 x 10 39 erg s^ 1 found in and around elliptical galaxies were 
consistent with background expectations, with a couple of 
notable exceptions. Hence this detection is quite remark- 
able, particularly at a luminosity of ~ 3 x 10 41 erg s _1 which 
makes it the second most luminous ULX can didate known 
after ESO 243-49 HLX-1 JFarrell et alj|2009h . We caution 
that there is a non- negligible chance of it being an unidenti- 
fied background object. However, we note that its location 
- close to a peculiar SO galaxy that possesses a dust lane, 
somewhat like Cen A but with no indication of nuclear activ- 
ity - is also interesting, particular ly as ESQ 243-49 H LX-1 
shares a very similar environment IjFarrell et al.ll2012i ). This 
object therefore deserves further study to elucidate its na- 
ture. 



4.2 X-ray characteristics 



In Section 13.21 we investigated the X-ray spectra of the ob- 
jects in this sample, finding a general preference for a power- 
law-like spectrum, which was intrinsically hard. We place 
these results in the context of other work on ULXs in Fig. [6] 
where we compare the photon indexes of the power-law fits 
to the objects in our sample to three more samples of ULXs 
fitted with power-law continua, all as a fun ction of luminos- 
ity. Th e three comparator samples are fromi lGladstone et al.l 
l|2009l ). containing 12 of the best quality ULX datasets in the 
XMM-Newton archive, and covering the 10 39 - 10 40 erg s _1 
luminosity r egime; the larg e r sam ple of Chandra datasets 
analysed by iBerghea et al.l (J2008), ranging from 10 39 to 
5 x 10 40 erg s _1 in luminosity; and the large sample of ULX s 
with XMM-Newton data presented by IWinter et al.1 (2006); 
although after removing their correction for absorption, we 
find the observed luminosities of many objects in the Winter 
et al. sample to lie significantly below 10 39 erg s _1 . 

Fig. [5] suggests a hardening of the average photon in- 
dex with increased observed luminosity; we investigate this 
further in Table [8] where we calculate the weighted mean 
value for the photon index for various samples shown in 
Fig. [6] This clearly shows that the average spectrum of 
ULXs is harder at high e r lum inosities, as originally sug- 
gested by IBerghea et al.l (|2008r i. Indeed, the effect appears 
quite pronounced; the two samples comp osed of sources be- 
low 10 40 erg s -1 ([Gladstone et ah 2009, and the less lumi- 
nous IBerghea et al.ll2008l objects) have weighted mean pho- 
ton indexes (r) > 2, whereas the more luminous samples 
are significantly harder with (r) ~ 1.6. This is clear ev- 
idence of a difference in the properties of the most lumi- 
nous ULXs from the more ordinary ULX population. Inter- 
estingly, similarly hard spectral indexes have been seen in 
two well-studied HL Xs: M82 X-l with e.g. V = 1.55 ± 0.05 
jKaaret et alj|2009l). and Cartwheel N 10 with Y = 1.75 ± 
0.25 |Wolter. Trinchieri fc Colpill2006l ). However the most 
luminous ULX, ESO 243-49 HLX-1, appears somewhat dif- 
ferent; in its bri ghtest phases it i s distinctly soft, with 
e.g. T = 3.4 ±0.2 (|Godet et aH^OQgT FH althou gh it hardens 
significantly to T ~ 2 when at fainter fluxes (jFarrell et al.1 



14 The best quality X-ray spectral data currently available for 
ESO 243-49 HLX-1 at its peak luminosities show its spectrum to 
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Figure 6. The relationship between power-law photon index V and observed luminosity for ULXs. Data points from our sample are plotted 
only where sufficient data was available for spectral fitting, and are plotted against the observed 0.3-10 keV luminosity, calculated using 
the power-law spectral fit and assuming that sources are at the distance of the identified host galaxy. All error bars are la uncertainty 
regions. Different sources are identified by the symbols shown in the key (top ri ght). Three samples o f lower luminosity ULXs are included 
for comparison (smaller data points). Data are taken from: purple plus sign — iBergh ea et al .l J2008I) . where luminosities were calculated 
from spectral parameters using a dummy respon se in XSPEC; blue fille d circles - Gladstone et al] J20091). luminosities were obtained 
from a coupled disc-corona model (DKBBFTH. iDone fc Kubotall2006i) ; red crosses - IWinter et al.l (2006), unabsorbed luminosities 
were converted to observed luminosities using PIMMS and assuming the power-law spectral fit. Typical error bars for each of the three 
comparison samples are shown in the same order from left to right in the top left corner. 



Table 8. Mean photon index for various ULX samples 



Sample 



U'} 



Glads tone et al. 2.19 ±0.01 

Berghea et al.. < 10 40 erg s" 1 2.10 ± 0.07 

Berghea et al.. > 10 40 erg s" 1 1.54 ± 0.06 

This sample (all) 1.72 ± 0.07 

This sample, > 5 X 10 40 erg s" 1 1.61 ± 0.07 



l201ll ; IServillat et alj I20T1I ) , such that it appears similar to 
the sources presented here (see section 15. ip . 

We also previously examined evidence for intra- 
observation variability in our objects. Here, we compare 
the results to well-constrained variability measurements 
from less-luminous ULXs. Specifically, in Fig. [7] we com- 
pare the fractional variability over the 0.3-10 keV range 
for our sample of luminous ULX candidates to the frac- 
tional variability of a selection of lower luminosity ULXs, 



be dominated by a c ool (fcT; n < 0.26 keV) disc component, see 
IServillat et alj <201ll) . 



on time scales of both 200 s and 2 ks. The comparison 
sample were selec ted for having b oth a well constrained 
power s pectrum in Heil et al. ([20091 ) . and a well-constrained 
flux by Idadstone et al] 1 20091 ) ■ Their fractional variabili- 
ties were calculated by integrating the power-law an d bro- 
ken power-law power spectra from iHeil et al. (120091) . and 
norma lising by flux using the values from Idadstone et al] 
(2009). As Fig. [7] demonstrates, the variability levels of 
the less luminous ULXs are low (typically 2 per cent frac- 
tional variability or less where well constrained), consis- 
tent with many r eports of low shor t -term variability lev- 
els in ULXs (e. g. ISwartz et afll2004 iFeng fc KaareH I2OO6I ; 
IHeil et al] 2009). In contrast, the measurements and con- 
straints on the higher luminosity objects are consistent with 
higher levels of variability, ~ 10 per cent or more fractional 
variability in most cases, although most error bars are rela- 
tively large due to the poorer data quality in these objects. 
This, then, tentatively points to a second observational dif- 
ference between 'normal' ULXs and their higher luminosity 
cousins; more short-term variability in the more luminous 
objects. 

In Fig. [8] we show the changes in the luminosity of the 
objects with time. It is immediately obvious the data is ex- 
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Figure 7. Fractional variability amplitude (F va r) shown as a function of 0.3—10 keV luminosity. This is plotted for time resolutions of 
200 s (top) and 2 ks (bottom). The upper limits shown are 2<r significance levels in all cases, and error bars are la confidence regions. Each 
object from the luminous ULX sample is labelled by a different symbol; a key is provided in the top left of each panel. Lower luminosity 
ULXs (small crosses) are shown for comp arison, these data points were obtai n ed by the integration of fits to the power spectra for a 
selection of sources from iHeil et al.l (120091 ') . Count rates from Idadstone et al.l ( 120091 ') for the same observations were used to calculate 
F var . The comparison sample includes (from left to right) - NGC 2403 X-l, NGC 1313 X-l, NGC 1313 X-2, NGC 5204 X-l, Ho IX X-l, 
NGC 4559 X-l and Ho II X-l. 
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Figure 8. Long term light curves of the objects in our sample. We show the observed 0.3-10 keV luminosity, plotted against observation 
date. Filled symbols correspond to observations with greater than 100 counts, where the fitted absorbed power-law models were used 
to calculate luminosity. Unfilled symbols correspond to observations with insufficient counts for spectral fitting, where instead spectral 
parameters were assumed as per Section 3.4. The low luminosity data point for Src. 3 corresponds to the most luminous of the resolved 
sources (Src. 3a) in the Chandra observation. 



tremely sparse, with no more than 2-6 measurements per 
object over a baseline of up to a decade; hence any con- 
clusions made on their long-term behaviour must be very 
tentative as we simply cannot tell how the objects behaved 
when they were not being observed. With this caveat in 
mind, we can tentatively sort the objects into two groups: 
firstly three objects that appear to fade substantially from 
an initial peak in luminosity, down to more 'usual' ULX 
luminosities (< 2 x 10 40 erg s" 1 ). This includes the only 
HLX candidate with a spiral host, Src. 1 in NGC 470, the 
XMM-Newton detection subsequently resolved by Chandra, 
Src. 3, and Src. 10 (which we discuss further below). The re- 
maining five objects all sustain their extreme luminosities, 
within a factor ~ 2, most notably Src. 7 which is detected at 
a luminosity ~3x 10 41 erg s _1 in both epochs. Given the 
sparse data available from our sample, and the lack of exten- 
sive monitoring of less luminous ULXs, we are also able to 
do little to compare the long-term variability properties of 
our objects to other ULXs. Recently Swift monitoring pro- 
grammes have been undertaken for a small handful of ob- 
jects, that appear persistently bright with typical flux vari- 
ations of factor s up to ~ 5 over the course of the moni toring 
campaigns fe.g. lKaaret fc Fendl2009l ; lGrise et al]|201Ch . The 
HLX candidate M82 X-l also appears to vary similarly in 
RXTE monitoring, altho ugh it also displays a periodic sig- 
nal on 62 day timescales (|Kaaret fc Fengll2007l ). The bright- 
est HLX, ESO 243-49 HLX-1, again appears different - its 



variability resembles an outburst cycle in an ordinary X-ray 
binary, with a fast rise to peak luminosities followed by a 
slow decay in flux, and a 375 ± 13 day gap between two 
outburst peaks l|Lasota et al.ll201ll ; IServillat et akll201ll ). 



Further inspection of Fig. [5] shows that several of 
the objects display significant X-ray spectral variability, 
in this case represented by a change in photon index at 
different flux levels (or, indeed, at different observation 
epochs with similar flux levels, as is the case for Src. 
8). However, these changes are not consistent across the 
small sample, with some objects hardening with increased 
flux (e.g. Src. 2), and others softening (e.g. Src. 1). This 
at least is consistent with the varied picture seen across 
the population of lower luminosity ULXs, where many 
different behaviours are observed with changing luminos- 
ity 



Feng fc Kaaret 2009; Kajava & Poutanenl [20091; 



Grise et all l20ld ; IVierdavanti et all l20ld ; IServillat et all 
20111 ). Perhaps more revealing are the situations in which 
the preferred spectral model changes with luminosity, as 
they may indicate direct analogies with known spectral state 
transitions. We discuss this in the following section. 
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5 DISCUSSION: ON THE NATURE OF THE 
MOST EXTREME ULXS 

This work has highlighted that there are possible differences 
between the standard population of lower luminosity ULXs, 
and those most luminous objects that we study here. Most 
notably, the brightest objects are spectrally harder, well de- 
scribed by power-law continua with photon indexes F ~ 1.7, 
and there is tentative evidence they display higher ampli- 
tude flux variability, with typical F var of ~ 10 per cent on 
timescales of 200 - 2000 seconds. But what does this mean 
physically for the objects in our sample? 



5.1 Intermediate-mass black holes? 

Perhaps the most interesting question we can ask is: 
do these objects finally offer evidence for a population 
of IMBHs underlying some ULXs? As discussed earlier, 
most super-Eddington models for stellar-mass black holes 
strug gle to explain the high luminosities of this sample 
(e.g. IZampieri fc Roberts! 120091 ). and there are suggestions 
that a different class of objects is required to e xplain sources 
abov e the X-ray luminosity function break (I Swart z et al.l 
I2011T ). It is therefore quite a straightforward step to infer 
that we may be looking at IMBHs. 

Crucially, this inference does appear to be supported 
by the observational characteristics of the sample. In par- 
ticular, the X-ray spectra and timing properties are consis- 
tent with the hard (forme r ly low /hard) state as described by 
iMcClintock fc Remillardl (|2006f ). sharing both a power-law 
spectrum (with T ~ 1.7) and a fractional variability level 
consistent with ~ 10 per cent. Given that the hard state 
is typical only found at Eddington fractions of 10 per cent 
or less, this immediately constrains the masses of our black 
holes to be in the IMBH regime. We calculate individual 
lower limits for each object, based on the highest observed 
luminosity with a power-law-like spectrum and the 10 per 
cent Z/Edd limit, in Table [5] as mass Mi. This implies black 
hole mass lower limits of 1.9— 27x 10 3 Mq across our sample. 

One issue with the timing result is that 
IMcClintock fc Remillardl J2006h note that 10 - 30 per 
cent fractional variability is expected in the 0.1 - 10 Hz 
range in the hard state, whereas we only probe frequencies 
as high as ~ 0.1 Hz in the best of our data. However, if we 
are dealing with large black holes we would expect their 
variability timescales to increase proportionally to thei r 
mass for the same accretion rate (jMcHardv et al.l 12006): 
hence a 1000 Mq black hole would vary on 100 times longer 
timescales than a 10 Mq black hole. So, the 0.1 - 10 Hz 
section of the PSD for a stellar mass black hole scales to the 
10~ 3 — 0.1 Hz regime, that we probe here, for a ~ 1000 Mq 
IMBH. 

An important test of whether these objects really are 
in the hard state would be to observe them transit into 
other states that are also seen at sub-Eddington luminosi- 
ties in Galactic binaries. Such t ransitions have been sug- 
gested for ESQ 243-49 HLX-1 l|Servillat et all 1201 ll ) and 
M82 X-l (|Feng fc Kaaretl |201Ch . with in both cases the 
high luminosities at which these transitions occur provid- 
ing a strong argument for IMBHs. Interestingly, ESO 243-49 
HLX-1 appears to possess a hard power-law X-ray spectrum 
at similar luminosities to the objects in this sample. We po- 



tentially see a similar transition in one case: Src. 1 is far 
better described by a MCD, representative of the thermal 
dominant state, than a power-law at its highest luminosity 
(~ 10 41 erg s _1 ), and vice versa in subsequent lower lumi- 
nosity datasets (< 2.5 x 10 40 erg s _1 ). Unfortunately, the 
timing data is of insufficient quality to confirm these states 
(with the little or no variability expected in the thermal 
dominant state potentially providing a distinction). How- 
ever, the temperature of the disc appears highly anomalous 
for an IMBH, at ~ 1 keV, making the identification of a 
classical transition for an IMBH very dubious in this case. 

One prediction that we can make is that if these are 
IMBHs in the hard state, then they should be emitting a 
steady radio jet. If so, the peak observed X-ray luminos- 
ity can be used in conjunction with the radio luminosity to 
provide limits on the black ho le mass via the so-called X- 
ray-r adio fundamental plane (|Merloni. Heinz fe di Matted 
|2003J). We have attempted to do this using the cur- 
rently available data. To this end, the VLA FIRST 
ilBecker. White fc Helfandlll995h and NVSS l|Condon et all 
ll998T ) 1.4 GHz radio surveys were searched at the source 
positions for any radio counterparts. However, none were 
detected, so the radio catalogue completeness limits were 
used as an upper-limit to the radio flux instead. A flat radio 
spectrum was assumed to extrapolate the radio flux from the 
observed band to 5 GHz, and 2-10 keV X-ray fluxes were 
extracted from the highest 0.3-10 keV flux of each source 
using pimms and spectral parameters from the appropriate 
power-law fit. Before presenting the limits we highlight two 
important caveats: firstly that the fundamental plane is not 
calibrated for IMBHs/ULXs, which may lead to some un- 
certainty; and secondly we are severely limited in what con- 
clusions we can draw by the lack of simultaneous X-ray and 
radio observations. 

Nevertheless, we present the upper limits on the IMBH 
masses in Table [9] as M2 . For a number of the sources the 
fundamental plane estimated limits on black hole mass are 
several orders of magnitude greater than the hard state mass 
limits, and in the SMBH mass regime. For these sources 
large improvements in radio detection limits would be nec- 
essary to use the X-ray-radio fundamental plane to test for 
IMBHs. However, for the closest source (Src. 8) the upper 
limit is already in the IMBH regime, and only ~ an order of 
magnitude above the lower limit from the hard state identi- 
fication. For two more objects (Src. 2 and Src. 5) a modest 
improvement of less than an order of magnitude would begin 
to constrain the mass within the IMBH regime. 

However, there are problems with interpreting these ob- 
jects as IMBHs. One such issue is the location of the ob- 
jects, the majority of which are found in the same star form- 
ing environments as lower-luminosity ULXs. While the ne- 
cessity to form only one such object in any galaxy neatly 
sidesteps the strong arguments for the bulk of the pop- 
ulations of ULXs in these regions ha ving to be predom i- 
nantly stellar-mass objects (|Kindl200 J: iMapelli et al.[ [2008), 
one still has to both form these very massive objects, 
and then capture a high mass star, in order to appear 
as a ULX. Even so, processes such as the assembly and 
subsequent collapse of a very massive star at the centre 
of a young, massive star cluster have been suggested as 
viab le IMBH formation channe ls in star forming regions 
(e.g. iPortegies Zwart et all 120041 ). and may provide a natu- 
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Table 9. Limits on the black hole mass for an IMBH from hard 
state Eddington arguments and the X-ray - radio fundamental 
plane. 



Source ID 


M!/M a 


M 2 /M 6 


Src. 1 


> 1900 


< 


2 x 10 6 


Src. 2 


> 3500 


< 


2 x 10 5 


Src. 3 


> 4700 


< 


2 x 10 6 


Src. 5 


> 6900 


< 


5 x 10 5 


Src. 7 


> 2.7 x 10 4 


< 


2 x 10 8 


Src. 8 


> 3200 


< 


3 x 10 4 


Src. 9 


> 5200 


< 


3 x 10 7 


Src. 10 


> 4700 


< 


4 x 10 6 



Notes: "Black hole mass limit from Eddington arguments, assum- 
ing that the source is in the hard state and accreting at < 10 per 
cent of the Eddington rate when at maximum luminosity with a 
power-law spectrum. For seven of eight objects this is the highest 
luminosity seen; for Src. 1 its highest luminosity spectrum is far 
better described by a MCD, so we use the second highest lumi- 
nosity dataset for this object. b Black hole mass derived from the 
X-ray-radio fundamental plane. 



ral environment in which to capture a star to f eed the IMBH 
IjHopman. Portegies Zwart fc A lexander 2004). However, in 
this case the lack of particularly bright and massive stellar 
clusters coincident with the ULXs may be an acute problem, 
even more so than the l ack of such coinc idences noted for 
less luminous ULXs fe.g. lZezas et alj|2002r ). Alternatively, if 
the IMBHs are not newly formed, the small number of such 
objects we observe may be commensurate with the predicted 
difficult ies of capturing a young st ar into a close orbit by an 
IMBH (JMadhusudhan et al.l 120061 ). presuming such objects 
could find their way into a star forming region. 

For the one object not readily related to star formation, 
Src. 7 in the SO galaxy IC 4320, assuming the identifica- 
tion as a hard state binary is correct then an older IMBH 
must be required. In this case the possibi lities include a 
prim ordial Population III star remnant (e.g. iMadau fc Reesl 
l2001f). or the stripped nuc l eus of an accreted dw arf galaxy 
fe.g. lKing fc Dehnen|[2005l ; iBellovarv et al.ll2010l ), although 
the requirement to accrete from a massive stellar reservoir 
likely remains. Again, we note the only established HLX 
with a luminosity in excess of Sr c. 7, ESQ 243-49 H LX-1, is 
also associated with an SO host (JFarrell et al.ll2012l ). 



5.2 The most extreme super-Eddington objects? 

While there are possible means of explaining the most lumi- 
nous ULXs in star forming regions as IMBHs, as discussed 
above, it is incontrovertible that this is the same environ- 
ment as that in which the majority of lower luminosity ULXs 
are found. We should therefore explore whether it is possible 
that these most luminous objects are driven by a more ex- 
treme version of the super-Eddington processes now thought 
to power the majority of the lower luminosity ULXs. 

There are possible scenarios in which we might explain 
most of these luminous objects by super-Edd i ngton emis- 
sion. Firstly, we note that lZampieri fc Roberts! ([2009) argue 
that the black holes underlying some ULXs in low metal- 
licity regions may be formed from the collapse of a single 



low-metallicity star which, by virtue of reduced wind-driven 
mass loss during its life time, can leave a black hole rem- 
nant as large as ~ 90 Mq (see also e.g. iBelczvnski et al.l 
|2010| for mass transfer binaries). A combination of this high 
mass and similar maximal Eddington rates as are inferred in 
some stellar mass ULX models could raise the maximum lu- 
minosity for stellar remnant ULX models to ~ 10 41 erg s _1 . 
In support, we note that the apparent locations of many of 
the ULXs studied here, on the outskirts of their host galaxy 
(cf. Fig. [2j , may indeed be consistent with the lower metal- 
licity regions required to form such objects. Secondly. iKind 
(2009) suggests that the degree of beaming due to geometric 
collimation by the inflated inner disc and outflowing wind 
at high accretion rates, that constrains the escape path of 
the X-rays from the central region, scales directly propor- 
tionally to the accretion rate. Hence, extreme luminosities 
may be explained by a high beaming factor (which may also 
help explain the rarity of such objects, as significant num- 
bers could be present, but would not be visible outside their 
narrow beam). 

However, both these scenarios would still have to ex- 
plain why the X-ray characteristics of these ULXs ap- 
pear so similar to the very sub-Eddington hard state, and 
not akin to the ultral uminous state spectra described by 
Idadstone et al.l (|2009i ). One possibility is the low quality of 
the data in most cases masks both the soft excess and the 
subtler high energy break of the ultraluminous state, and 
the underlying spectra are similar to those observed in Ho 
IX X-l or NGC 1313 X-l bv ldadstone et all where the best 
fitting power-law continua have indexes ~ 1.6 — 1.85. How- 
ever, these observations both have particularly low fractional 
variability (cf. Fig. [TJl , unlike the objects in the higher lu- 
minosity sample, which appear variable on short timescales 
(albeit with large errors, or only upper limits in some cases) . 
One way around this would be if the higher luminosity ob- 
jects are more prone to displaying extrinsic variability. This 
is suggested by Middlct on et al.l ([2011) to explain the vari- 
ability of NGC 5408 X-l, where it could originate in blobs 
of outflowing material driven off the surface of the accre- 
tion disc by the extreme radiation pressures from its central 
regions, that cross the line-of-sight to those regions. Better 
data is required to resolve these issues; but it is notable that 
in the one dataset within this sample with reasonably good 
data quality (Src. 8 in NGC 5907) an ultraluminous state 
spectrum is observed. However, we note that the distance to 
NGC 5 907 is uncertain, and recent estimates (e.g. d — 13.4 
Mpc bv lTullv et al.l l2009) place the host galaxy near enough 
that the luminosity of Src. 8 is reduced by a further 20 per 
cent, such that it is consistent with the top end of the 'nor- 
mal' ULX population. It is therefore questionable whether 
it should be included within this sample. Even so, its hard 
power-law-like spectrum and detected fractional variability 
provide important testimony that an ultraluminous state ob- 
ject can possess similar observational characteristics to the 
sub-Eddington hard state in moderate quality data. 

One further mass limit could be placed assuming 
that the sample contains hi gh a ccretion rate objects. Both 
iGladstone fc Roberts! (|2009l ) and lYoshida et al.l feoich note 
that the low luminosity ULXs in NGC 4490 and (in the 
latter paper) M51 appear to transit between a power-law 
spectral state at ~ 10 39 erg s -1 to a MCD-like m odel at 
slightly higher luminosities. IGladstone fc Roberts! suggest 
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this could be a transition from a sub-Eddington steep power- 
law state toa~ Eddington rate state; if so the highest lumi- 
nosity in which a power-law spectrum is observed provides 
an estimate of the black hole mass from the requirement 
Lx/LEdd < 1. However, without seeing the transition lu- 
minosity to a MCD spectrum itself this measurement does 
not provide strong constraints (as with further observations 
we might see the power-law at higher luminosities); for the 
sources in the sample this estimated upper limit is in the 
range ~ 200-2700 Mq, with the only possible constrained 
limit (from Src. 1) as < 800 Mq from the observed luminos- 
ity of the observation where a MCD fit is much preferred. 

5.3 Contaminants? 

A final point to consider is whether any more of our 
objects could be background (or, plausibly, foreground) 
contaminants . A calculation based on that performed by 
IWalton et al.l (2011b) was used to establish the expected 
contamination of this sample by background objects; this 
found that three were expected out of the ten selected 
sources with luminosities above 5 x 10 40 erg s _1 within 100 
Mpc. In the analysis we were able to identify two back- 
ground QSOs near elliptical galaxies, including one new 
2 = 3.25 QSO near NGC 4874. The calculation predicts one 
more background object, likely to be associated with a spi- 
ral galaxy, but no such object has been identified at present. 
We note that we would be unlikely to identify such an object 
from its X-ray properties alone. In Table A.l of Appendix A 
we present X-ray spectral fits to the QSO datasets, using the 
same models as for the ULXs (cf. Table [SJ). Although Src. 
4 does appear somewhat softer than the luminous ULXs, 
errors on its spectral parameters are large. Large errors also 
dominate the measurements for Src. 6, and we were not 
able to derive short-term variability characteristics for either 
source due to the paucity of data. We would not, therefore, 
have been able to distinguish these objects from ULXs with- 
out the optical follow-up data. We therefore caution that 
without further optical follow-up (which, as the case of Src. 
6 shows, must go several magnitudes deeper than the DSS) 
we cannot rule out a QSO nature for the several ULXs cur- 
rently lacking follow-up data. 

Further possible contaminants include foreground ob- 
jects and other rare types of objects within the host galax- 
ies. We can at least rule out foreground stars as none are 
evident in the DSS images in FigO However, one possibility 
that should be mentioned is that of recent supernovae in the 
host galaxies, whic h are not accou nted for in the contamina- 
tion calculation of IWalton et al.l (2011b). Such supernovae 
have been seen with typical X-ray luminosities in the range 
of 10 37 -10 41 erg s _1 and fade over a variety of timescales, 
e.g. SN 1988Z decreased by ~ 5 over 2 years, whilst SN 
1980K decreased by a similar factor bet ween days 35 and 
82 after outburst Ijlmmler fc Lewin 2003). Interestingly, the 
X-ray spectra of some objects have been see n to soften as 
the fl ux decreases, for example in SN 1999em (|Poolev et al.l 
120021 ). where the early X-rays are thought to originate in the 
circumstellar gas shocked by the initial blast wave, whereas 
the later soft emission co mes from the lower temperature 
rever se shock region (see iFransson. Lundqvist fc Chevalierl 
Il996l for details). We do see several objects in the sample 
whose X-ray emission appears to fade in the years following 



their initial detection, including one object (Src. 10) that 
also softens dramatically. However, each of these objects 
shows some degree of fractional variability within at least 
one observation, consistent with a compact object and re- 
jecting the recent supernova hypothesis. This, and the lack of 
any identified optical supernovae associated with the ULXs 
in this sample, argues that we do not have any such objects 
amongst our sources. 

A second rare class of objects that may be associated 
with the host galaxies are recoiling SMBHs, as has been 
suggested for the H LX candidate CXO J122518. 6+144545 
IjJonker et alj|2010j ). Such objects are predicted to form via 
the binary coalescence of SMBHs during galaxy mergers, 
where anisotropic gravitational wave emission imparts a kick 
of u p to several tho usand km s _1 on the new, single SMBH 
(see lKomossall2012i and references therein). This kick is suf- 
ficient to displace the SMBH from the centre of its host 
galaxy, and in some cases even expel it completely. Such 
objects might appear both X-ray and optically luminous, as 
they are predicted to retain a reservoir of gas and stars to ac- 
crete from; a key indicator of their SMBH nature would then 
be broad emission lines in their optical spectra. However, we 
again regard the chances of detecting such an object in our 
sample to be low; none of our spiral galaxies are undergoing 
(or appear to have recently undergone) the major merger 
necessary to create such an object. Interestingly, the only 
early-type galaxy in our sample, IC 4320, possesses a dust 
lane, a characteristic tha t has been linked to minor galaxy 
mergers in such systems (|Shabala et al.ll201lr ). It also hosts 
the most luminous ULX in our sample. The detection and 
characterisation of the optical counterpart to this object is 
therefore of great interest. 



6 CONCLUSIONS 

Even in the absence of definitive black hole mass measure- 
ments for ULXs, various arguments have resulted in the wide 
acknowledgement that the majority of these extraordinary 
objects are most probably powered by super-Eddington ac- 
cretion onto stellar remnant black holes. However, the lu- 
minosities that can be produced in this manner are limited 
by both a maximum Eddington ratio and maximum stellar 
remnant black hole mass, which complicates the unification 
of the brightest ULXs with the less luminous population. 
These particularly luminous, extreme ULXs may instead re- 
main as potential IMBHs, with their rarity consistent with 
arguments as to the difficulties in producing stable high ac- 
cretion rates in many such objects. 

Here we have presented a study of 10 candidate extreme 
ULXs. We were able to identify 2 as background QSOs, re- 
sulting in a sample of 8 good candidates for further study. 
By examining the X-ray emission characteristics of this sam- 
ple we note that these extreme sources do indeed appear to 
be distinctly different from the lower luminosity ULXs. Not 
only are they typically spectrally harder, they are also con- 
sistent with being more variable (albeit at low statistical 
significance in many objects). This combination of spectral 
and timing characteristics is highly suggestive of the sub- 
Eddington hard state, which given the luminosities involved 
would strongly argue for large (10 3 -10 4 M ) IMBHs. This is 
however not entirely conclusive, as accretion at the highest 
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plausible super-Eddington ratios onto the most massive stel- 
lar remnant black holes could potentially produce ULXs up 
to ~ 10 erg s _1 . Further high quality X-ray observations 
are therefore required to test this, with the capability to 
detect ultraluminous state-like signatures, especially for the 
sources that lie in the transition zone between the brightest 
local ULXs (~3x 10 40 erg s" 1 ) and the HLX luminosity 
range. 

The sample of 8 candidate ULXs includes 2 candidate 
HLXs. One is a transient HLX, that subsequently fades to 
high ULX luminosities; but the other appears persistent be- 
tween two snapshots separated by 8 years, at a luminos- 
ity that is second only to ESO 243-49 HLX-1 amongst the 
known HLXs. These join the very small number of HLXs 
identified in the literature. As they are particularly difficult 
to explain as stellar remnants, they must be regarded as 
some of the best known IMBH candidates. Given the cos- 
mological and astrophysical importance of IMBHs, HLXs 
are highly compelling targets for further observations with 
current and future missions, to both better constrain the 
presence of a massive black hole and better understand the 
astrophysics of these remarkable objects. 
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Figure Al. Gemini GMOS-N longslit spectrum of the optical counterpart to Src. 6, taken using the R150 grating and the 0.75 arcsecond 
slit. It was suggested that this source may be a bright ULX in a UCD in the Coma cluster, however lines identified in this spectrum 
(shown by dashed red lines; see text for details) instead identify it as a QSO at redshift £S 3.25. 
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d 0.3— 10 keV luminosity (xlO 40 erg s — 1 ), estimated from a power-law spectral fit, and assuming the object is at the distance of the 
foreground galaxy. e 0.3-10 keV luminosity (xlO 44 ergs -1 ), estimated from a power-law spectral fit, and assuming the object is at 
the distance derived from its optical spectral redshift measurement. -^This observation had insufficient counts for spectral fitting, so a 
power-law spectrum with the same spectral index and absorption column as observation 0112271001 was assumed when estimating the 
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APPENDIX A: THE TWO IDENTIFIED 
BACKGROUND CONTAMINANTS 

In the course of our work we were able to identify two ex- 
treme ULX candidates as background contaminants. These 
were Src. 4, initially associated with NGC 4065 by Walton et 
al. (2011) and identified as a background QSO on the basis 
of SDSS spectroscopy; and Src. 6, initially associated with 
NGC 4874 but later identified as a background QSO using 
Gemini GMOS spectroscopy. In Table [ATI we present the re- 
sults of the simple X-ray spectral fits to the data for these 
objects, for comparison with the presumed bona fide extreme 
ULXs analysed in the body of the paper (see Section 5.3 for 
discussion related to this comparison). Here we also present 
the Gemini GMOS spectrum of the optical counterpart to 
Src. 6 in Figure IA1I The broad emission lines present at 
~ 5177, 5265 & 8076 A are identifiable with the Lyman a 
1216 A, NV 1240 A & CIII 1909 A lines (respectively) at a 
redshift of z « 3.25, demonstrating this is very likely to be 
a high redshift QSO. 

This paper has been typeset from a TjtX/ P/TpX file prepared 
by the author. 
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